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Nomenclature 


a, b, c et d : coefficients semi-empiriques (polarisation de concentration) 

A k : porosité de la membrane 

C; (x) : concentration de l’ion (ou du soluté non chargé) i dans le couche de 

polarisation 

C; (x) : concentration de l’ion (ou du soluté non chargé) i à l’intérieur des pores 

c| nt : concentration de l’ion (ou du soluté non chargé) i à l’extérieur de la membrane 

et à l’interface membrane / solution externe 

c‘ nt : concentration de l’ion i à l’intérieur de la membrane et à l’interface membrane 

/ solution externe 


c : 5 : concentration d’alimentation du soluté non chargé i (polarisation de 

concentration) 

DCE double couche électrique 

D,: coefficient de diffusion à dilution infini de l’ion (ou du soluté non chargé) i 

d h : diamètre hydraulique de la section d’écoulement 

D se i: coefficient de diffusion effectif du sel 

e: charge élémentaire 

F: constante de Faraday 

int: peut être remplacé par 0 ou Àx selon l’interface considérée 

I : force ionique de la solution 

Io, Ii : fonctions de Bessel modifiées 
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Nomenclature 


jp flux molaire de l’ion 

J s : flux de soluté transféré à travers la membrane 

J v : flux volumique du perméat 

J w : flux volumique de l’eau 

k: coefficient de transfert de masse dans la couche de polarisation 

k’ vecteur d’onde 

k|j: constante de Boltzmann 

Ko, Kp fonctions de Bessel modifiées 

K IX : coefficient hydrodynamique rendant compte de l’effet de la paroi du pore sur 

le transport convectif 

Ki.o: coefficient hydrodynamique rendant compte de l’effet de la paroi du pore sur 

le transport diffusif 

Kj s K; t : fonction hydrodynamique 

L : longueur de la membrane 

Lp : perméabilité hydraulique 

P: périmètre mouillé 

Pe: nombre de Péclet 

PHE : Plan d’Helmholtz Externe 

PHI : Plan d’Helmholtz Interne 

r; : rayon ionique de l’ion i 

q cav : rayon de cavité de l’ion i 

r s : rayon de Stokes de l’ion i 

r p : rayon du pore moyen ou demi-hauteur moyenne suivant la géométrie 
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Nomenclature 


R: 

constante des gaz parfaits 

Re: 

nombre de Reynolds 

Rfint : 

taux de rejet intrinsèque de l’ion (ou du soluté non chargé) i 

Rfobs • 

taux de rejet observé de l’ion (ou du soluté non chargé) i 

Rm : 

résistance hydraulique (intrinsèque) de la membrane 

S : 

aire de la section d’écoulement 

S 1/2 : 

sélectivité d’une membrane pour un composé 1 par rapport à un composé 2 

Sc: 

nombre de Schmidt 

Sh: 

nombre de Sherwood 

T: 

Ti: 

température 

taux de transmission de l’espèce i 

v : 

vitesse de circulation de la solution d’alimentation 

V: 

vitesse du fluide dans les pores 

Wi : 

énergie libre d’interaction incluant toutes les interactions de l’ion avec le 
milieu (énergie de solvatation) 

x: 

coordonnée 

X: 

densité de charge volumique effective 

Zi-' 

charge de l’ion i 


Symboles 

ô: 

épaisseur de la couche de polarisation 

Ax: 

épaisseur de la membrane 

AP: 

différence de la pression hydrostatique 
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Nomenclature 


AW' ,int • 

excès de l’énergie de solvatation normalisée due à l’effet de Born de l’ion i 
normalisé par k B T 

AW.’ int : 

i, image 

excès de l’énergie de solvatation normalisée due aux forces images normalisé 
par k B T 

A^cell : 

potentiel de cellule 

^diff : 

potentiel de diffusion 

A'Pq 1 : 

potentiel de Donnan 

AU/ 

^ * élect * 

potentiel de concentration 

A'P • 

potentiel de membrane 

£o: 

permittivité du vide 

e b : 

constante diélectrique de la solution à l’extérieur des pores 

£m- 

constante diélectrique du matériaux membranaire que constitue la couche 
active 

£ p : 

constante diélectrique des pores 

cf)| nt : 

coefficient de partage stérique de l’ion i 

y | nt . 

coefficient d’activité ionique à l’extérieur de la membrane 

— int 

Y, : 

coefficient d’activité ionique à l’intérieur de la membrane 

rp 

viscosité dynamique du fluide 

K 1 : 

longueur de Debye 

k j nt : 

h: 

coefficient de partage stérique / diélectrique de l’ion i 
rapport entre le rayon de Stokes de l’ion i et le rayon du pore 
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Nomenclature 


p. : potentiel électrochimique d’un ion i à l’extérieur des pores 

p. : potentiel électrochimique d’un ion i à l’intérieur des pores 

: potentiel chimique standard de l’ion i 

p" 11 : longueur d’écrantage effective réciproque normalisée 

p : masse volumique 

ao : densité de charge de surface 

op : densité de charge de la couche compacte 

ad : densité de charge de la couche diffuse 

'P: potentiel électrostatique local 

'Po : potentiel de surface 

*Pp : potentiel au PHI 

V P (I : potentiel au PHE 

'P 1111 : potentiel électrostatique local à l’extérieur des pores et à l’interface membrane 

/ solution externe 

^p "' 1 : potentiel électrostatique local à l’intérieur des pores et à l’interface membrane 

/ solution externe 
Ç : potentiel zêta 
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Introduction 


Depuis les quatre dernières décennies, les techniques de séparation mettant en œuvre 
des membranes suscitent d’importants efforts de recherche et de développement. Ces 
procédés sont de plus en plus utilisés dans des secteurs d’activités aussi nombreux que variés 
comme les industries chimique, agro-alimentaire, pharmaceutique, la biotechnologie, la 
production d’eau potable, le traitement des effluents industriels... Ce développement est lié 
aux avantages que les procédés à membranes présentent par rapport aux techniques de 
séparation conventionnelles. En effet, dans la plupart des cas : 

• l’opération de séparation est effectuée à température ambiante. Ainsi, un milieu 
thermosensible peut être traité sans qu’aucun de ses constituants ne soit altéré. Ceci est 
particulièrement important dans les industries agro-alimentaire et la biotechnologie. 

• le fractionnement se fait sans changement de phase, ce qui est très économique sur le 
plan énergétique. 

• le procédé est physique, donc l’installation est modulable et peut fonctionner en 
continu. 


Le développement des techniques membranaires est également lié à la sévérité 
croissante des lois environnementales qui a incité les industriels à investir dans des techniques 
de retraitement d’eaux usées. En matière de rejet d’eaux industrielles usées, les nouvelles 
directives européennes sont de plus en plus draconiennes : on parle de rejet zéro d’ici 2020 
(directives 2455/2001/CE et 2000/60/CE). D’autre part, dans les pays industrialisés, la 
consommation d’eau ne cesse de croître, ce qui génère des volumes d’eaux usées 
considérables, lesquels doivent être par la suite traités de manière à répondre aux exigences 
imposées par les règlements en vigueur. Le traitement des eaux usées, domestiques et 
industrielles, est devenu une priorité pour tout secteur industriel et toute communauté urbaine. 
Parmi les solutions possibles, les techniques membranaires occupent une place privilégiée 
grâce à leurs performances énergétiques et technico-économiques de plus en plus 
concurrentielles. 
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Introduction 


De façon générale une membrane peut être décrite comme une barrière sélective 
séparant deux milieux et permettant le passage préférentiel d’une espèce parmi d’autres sous 
l’effet d’une force de transfert. Les procédés membranaires permettent donc de concentrer, de 
purifier ou de fractionner une solution en contrôlant et en orientant les échanges de matière à 
travers la barrière sélective que constitue la membrane. 

Les procédés membranaires regroupent un grand nombre de techniques permettant de 
réaliser des séparations en phase liquide ou en phase gazeuse sous l’action de différentes 
forces de transfert (mécanique, électrique ou chimique). On y distingue les procédés 
baromembranaires, lesquels utilisent une différence de pression comme force motrice. Pour 
les séparations en phase liquide, ceux-ci sont au nombre de quatre : la micro-, l’ultra-, la 
nanofiltration, qui mettent enjeu des membranes possédant une structure poreuse, et l’osmose 
inverse qui utilise des matériaux denses. La microfiltration est utilisée pour séparer des 
espèces dont la taille est de l’ordre du micromètre (cellules animales, levures, bactéries...), 
l’ultrafiltration pour séparer des espèces de taille comprise entre ~2 et 100 nm (peptides, 
protéines, virus...) et la nanofiltration pour traiter des espèces de taille inférieure à 2 nm 
(sucres, colorants, sels....). Comme la nanofiltration, l’osmose inverse permet de filtrer des 
solutions contenant des espèces de très faible masse molaire. Elle est généralement utilisée 
lorsqu’une rétention totale des solutés est souhaitée. 

La nanofiltration est la plus récente des techniques membranaires. Son champ 
d’application se situe entre ceux de l’ ultrafiltration et de l’osmose inverse. Depuis sa récente 
implantation industrielle (années quatre-vingt), la nanofiltration a su trouver un nombre 
grandissant d’applications telles que le recyclage d’eau dans les effluents de teinturerie, la 
production d’eau potable, la déminéralisation du lactosérum, la décoloration des sirops de 
sucre, le traitement des eaux de blanchiment de la pâte à papier (contenant des dérivés 
organiques halogénés), le recyclage des eaux de lavage dans l’industrie des traitements de 
surface, le recyclage de solution alcalines de nettoyage, le recyclage d’eau à partir d’effluents 
de teinturerie, la concentration d’antibiotiques, l’élimination de pesticides et d’herbicides... 

Comme le suggèrent les quelques applications énoncées ci-dessus, la nanofiltration 
peut se révéler dans certains cas comme une technique concurrente de l’osmose inverse 
permettant certes des taux de rejet moins élevés mais présentant l’avantage d’utiliser des 
pressions transmembranaires plus faibles que celles utilisées en osmose inverse tout en 
produisant des flux de solvant supérieurs. 
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Introduction 


A la différence de l’osmose inverse, elle présente aussi une forte sélectivité entre les 
ions mono- et multivalents. Ceci provient du fait que les membranes de nanofiltration sont 
élaborées pour la plupart à partir de matériaux possédant des groupements fonctionnels 
ionisables en surface. Les membranes de nanofiltration possèdent donc une certaine charge 
fixe en solution. Parce que les pores de ces membranes ont des dimensions du même ordre de 
grandeur que les distances caractéristiques d’action des forces électrostatiques, leur sélectivité 
vis-à-vis de solutés chargés ne dépend plus seulement des effets stériques (liés aux tailles 
relatives des pores et des solutés) mais également des interactions d’origine électrostatique 
incluant les forces coulombiennes et très probablement des phénomènes d’exclusion 
diélectrique. Pour comprendre et prévoir la sélectivité d’une membrane de nanofiltration, il 
est donc indispensable d’une part, de déterminer ses caractéristiques structurales (taille de 
pore moyen, porosité, épaisseur...) et sa charge de surface et d’autre part, de disposer d’outils 
fiables de modélisation reliant les caractéristiques d’un nanofiltre à ses propriétés de transfert 
afin de pouvoir simuler le comportement de tout nanofiltre vis-à-vis d'un fluide à traiter. 

Un tel modèle de transport, désigné sous le vocable SEDE (Steric, Electric and 
Dielectric Exclusion), a été récemment développé au laboratoire. La pertinence de ce modèle 
réside sur la prise en compte simultanée de trois types d’interaction : l’exclusion stérique, 
l’exclusion de Donnan et les effets diélectriques incluant à la fois les effets de Bom 
(provenant de l’écart entre les constantes diélectriques de la solution externe et de celle dans 
les nanopores) et les forces images (induites par la différence entre la constante diélectrique 
du matériau membranaire et celle de la solution dans les nanopores). Ce modèle repose sur 
quatre paramètres d’ajustement qui sont le rayon de pore moyen, le rapport épaisseur sur 
porosité, la densité de charge volumique à l’intérieur des pores (nombre de moles de charge 
fixes par unité de volume de pore) et la constante diélectrique dans les nanopores. Parmi ces 
quatre paramètres, la détermination expérimentale des deux derniers est extrêmement 
difficile, voire impossible, ce qui rend difficile l’évaluation des qualités prédictives du modèle 
en terme de taux de rejet. 

Pour contourner ce problème, il a donc été décidé de tester la fiabilité du modèle 
SEDE en s’intéressant à sa cohérence, c’est-à-dire en analysant à la fois le taux de rejet d’ions 
et une autre grandeur caractéristique du transport, à savoir le potentiel de membrane. Cette 
grandeur représente la différence de potentiel électrique s’établissant de part et d’autre d’une 
membrane sous l’action d’une différence de concentration. Pour cela, le modèle SEDE a dû 
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Introduction 


être étendu à l’étude de cette nouvelle grandeur. Le modèle a tout d’abord été appliqué à 
l’étude théorique de l’influence de différents paramètres tels que le rayon de pore, la densité 
de charge volumique à l’intérieur des pores et la constante diélectrique de la solution à 
l’intérieur des pores sur le potentiel de membrane dans le cas de solutions binaires. Il est à 
noter que l’influence à la fois des effets stériques, électriques et diélectriques sur le potentiel 
de membrane n’avait jamais été abordée jusqu’à présent. Le modèle a ensuite été utilisé pour 
déterminer les densités de charge volumique et les constantes diélectriques des nanopores 
obtenues d’une part, à partir des taux de rejet expérimentaux de sels simples et d’autre part, à 
partir des potentiels de membrane mesurés avec ces mêmes sels. L’existence de paramètres 
permettant de décrire simultanément les deux types d’expériences (taux de rejet et potentiel de 
membrane) peut être considérée comme une preuve de la cohérence du modèle SEDE. 

Enfin, le modèle a été étendu à l’étude du transport d’un mélange ternaire sous l’action 
d’un gradient de concentration. On s’est alors attaché à étudier l’influence des effets stériques, 
électriques et diélectriques sur le potentiel de membrane. Il est à noter que le choix d’un 
mélange ternaire présente l’avantage de pouvoir estimer la constante diélectrique dans les 
nanopores à partir des seules mesures de potentiel de membrane. 


Le présent mémoire comprend trois chapitres. Le premier rappelle les principes de 
bases de la filtration membranaire et permet notamment de présenter les principales propriétés 
des membranes vis-à-vis d’un fluide à traiter. 

Le deuxième chapitre est dédié à la description des modèles de transport à travers des 
membranes de nanofiltration. Après avoir présenté un historique des modèles existant, les 
bases du modèle de transport SEDE sont rappelées puis ce dernier est appliqué à la 
modélisation du taux de rejet ionique et du potentiel de membrane. 

Le dernier chapitre se présente sous forme d’un recueil de publications présentant les 
principaux résultats issus de ce travail. 
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1 Généralités sur les 
membranes de filtration 


Une membrane peut être définie comme une barrière sélective ou semi-perméable de 
quelques centaines de nanomètres à quelques millimètres d’épaisseur séparant deux 
compartiments, et permettant le passage préférentiel d’une espèce parmi les autres sous 
l’action d’une force de transfert (figure 1). Dans le cas des procédés de filtration 
baromembranaire, la force de transfert nécessaire à la réalisation de la séparation est une 
différence de pression hydrostatique appliquée de part et d’autre de la membrane. 

Solution à traiter 

Force motrice 


= Q/» Q/ Éi/ 



Membrane 

Figure 1 : Représentation schématique d’une membrane. Schéma de principe de 
fonctionnement d’une membrane permsélective [Szymczyk et al., 2007]. 


Ce chapitre rassemble des généralités sur les membranes de filtration et leurs 
propriétés vis-à-vis d’un fluide. 
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1 Généralités sur les membranes de filtration 


1 . 1 Les différents types de 

membranes de filtration 


Il existe une grande diversité de membranes qui se différencient les unes des autres 
par leur nature, leur structure et leur géométrie. Les caractéristiques chimiques et structurales 
de ces matériaux leur confèrent des performances de séparation clairement identifiées les 
rendant ainsi aptes à réaliser des séparations spécifiques dans divers secteurs de l’industrie 
(traitement des eaux usées, des effluents industriels, des fluides alimentaires...). 

Nous détaillons ici les différents types de membranes utilisées dans les procédés 
industriels de filtration baromembranaire. 

1.1.1 Nature chimique du matériau 

Les membranes utilisées dans les procédés de filtration baromembranaire peuvent être 
de nature organique ou minérale (tableau 1). On distingue généralement : 

• les membranes de l ere génération à base d’acétate de cellulose qui sont les premières à 
avoir fait leur apparition sur le marché ; Ce sont des matériaux possédant de bonnes 
propriétés de perméabilité mais qui souffrent de conditions limites d’utilisation relativement 
contraignantes en raison de leur sensibilité à l’hydrolyse chimique et aux températures 
élevées ; 

• les membranes dites de 2 eme génération qui sont constituées de polymères de synthèse 
(poly sulfones, polyamides aromatiques, polycarbonates, polymères fluorés) ; Ces 
membranes, possèdent des résistances chimique et thermique supérieures aux membranes 
cellulosiques et ont permis d’augmenter significativement le nombre d’applications des 
procédés membranaires ; 

• les membranes inorganiques qui ne sont apparues que plus récemment sur le marché ; 
Constituées à base de matériaux céramiques (alumine, dioxyde de titane, dioxyde de 
zirconium), elles sont plus résistantes mécaniquement, chimiquement et thermiquement que 
leurs concurrentes organiques, ce qui justifie leur utilisation dans des applications mettant en 
jeu des milieux agressifs et/ou non aqueux ; De plus, elles peuvent être stérilisées à la vapeur 
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1.1 Les différents types de membranes de filtration 


sans altération de leur structure, ce qui les rend attractives dans certaines branches du secteur 
agroalimentaire comme l’industrie laitière par exemple. 


Matériau 

Température 
de travail (°C) 

Gamme 
de pH 

Avantage(s) 

Inconvénient(s) 

Membranes de l eie 
génération à base d’acétate 
de cellulose 

30-40 

3-8 

Les 1 ères sur le 
marché. 

Bonne perméabilité. 
Nettoyage facile. 

Sensibilité (pH extrêmes, 
température, oxydants). 

Membranes de 2 eme 
génération : 

• polyamide aromatique 

• polysulfone 

• polyacrylonitrile 

• polyfluorure de 
vinylidène 

60-80 

60-100 

60 

130-150 

2-11 

1-13 

1-10 

1-13 

Bon marché. 
Bonne résistance. 

Fragile au contact avec 
les agents chlorés 
(nettoyage). 

Problème aux fortes 
pressions. 

Membranes céramiques : 

• A1 2 0 3 a 

• Ti0 2 

• Zr0 2 

>900 

350 

400 

1-14 

1-14 

1-14 

Meilleure résistance 
chimique, thermique 
et mécanique. 
Stérilisation à la 
vapeur possible. 

Prix élevé. 


Tableau 1 : Quelques exemples de membranes couramment utilisées. [Fievet et al., 2009]. 


1.1.2 Structure 

1. 1.2.1 Membranes poreuses et membranes denses 


Une membrane poreuse présente une structure semblable à une éponge avec des pores 
continus pouvant être interconnectés entre eux. Un pore est défini comme un interstice entre 
les constituants solides de la membrane. Dans le cas des membranes organiques, la structure 
poreuse résulte des espaces entre les chaînes de polymère tandis que dans le cas des 
membranes inorganiques les pores sont constitués par les espaces inter- granulaires. Même 
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1 Généralités sur les membranes de filtration 


s’il existe quelques types de membranes ayant des pores droits et de géométrie bien définie 
(cas des membranes "track-etched" obtenues par bombardement ionique suivi d’une étape de 
révélation des traces ou des membranes obtenues par oxydation anodique de l’aluminium), la 
plupart des matériaux commercialisés sont caractérisés par un réseau de pores tortueux et 
interconnectés. 

Selon le type de membranes, le diamètre moyen des pores peut varier du nanomètre 
jusqu’à quelques dizaines de micromètres. L’IUPAC (International Union of Pure and 
Applied Chemistry) distingue trois catégories de pores selon leur taille et recommande 
d’adopter la terminologie suivante : 

• micropores, pour des diamètres inférieurs à 2 nm 

• mésopores, pour des diamètres compris entre 2 et 50 nm 

• macropores, pour des diamètres supérieurs à 50 nm 

Il faut noter que cette terminologie est assez rarement employée par les utilisateurs de 
procédés membranaires du fait de la confusion qu’elle génère avec les dénominations des 
différents procédés baromembranaires (tableau 2). A titre d’exemple, en adoptant la 
terminologie recommandée par l’IUPAC, la nanofiltration fait appel à des micropores tandis 
que la microfiltration utilise des macropores. La technique d’ultrafiltration recouvre quant à 
elle les deux catégories : mésopores et macropores. 



Microfiltration 

Ultrafiltration 

Nanofiltration 

Osmose 

Inverse 

Diamètre de pore 

10-0,1 (im 

100-2 nm 

2-0,5 nm 

Membrane 

dense 

Rétention 

Microorganismes 
et particules 

Macromolécules 

Ions multivalents et 
molécules 

Ions 

Perméation 

Solvants et 
espèces dissoutes 

Solvants, sels et 
petites 
molécules 

Solvants, ions 
monovalents et 
petites molécules 

Solvants 

Pression 

transmembranaire 

0,1 à 5 bar 

0,5 à 9 bar 

4 à 20 bar 

> 20 bar 


Tableau 2 : Classification des procédés baromembranaires en fonction de la taille des pores 
et des espèces retenues [Les Cahiers du CFM n°4]. 
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1.1 Les différents types de membranes de filtration 


La porosité d’une membrane est définie comme le rapport du volume des espaces 
vides sur le volume total de la matrice. Cette grandeur dépend de la taille des constituants 
élémentaires de la matrice mais également du procédé d’élaboration de la membrane. Une 
membrane poreuse doit posséder une très bonne résistance mécanique tout en permettant un 
débit de perméation élevé. La séparation dépend des effets stériques (c’est-à-dire d’une 
discrimination en fonction de la taille de l’espèce à retenir), mais aussi des effets de charges 
électriques découlant de la nature chimique du matériau membranaire. 


Lorsque les interstices entre les constituants de la matrice membranaire ont une taille 
voisine de celle d’un ion solvaté, on ne parle plus de pores et on qualifie la membrane de 
dense. Elle est soit à base de polymère de haute densité se présentant sous forme de feuilles 
souples, soit à base de matériaux céramiques. Le transfert de matière à travers une membrane 
dense est régi par un mécanisme de solubilité des solutés dans la membrane et leur diffusion. 


1. 1.2.2 Membranes symétriques et membranes 
asymétriques 


Les membranes peuvent également être classées en fonction de la structure des 
différents matériaux qui la composent. On distingue : 

• les membranes symétriques (ou isotropes) dont la structure est homogène sur toute 
l’épaisseur (figure 2a). Etant donné que la résistance intrinsèque d’une membrane est 
proportionnelle à son épaisseur, il est donc nécessaire qu’une membrane symétrique soit la 
plus mince possible. 

• les membranes asymétriques (ou anisotropes) dont la structure varie d’une couche à 
l’autre (figure 2b). Celles-ci présentent une structure hétérogène. Il peut s'agir de membranes 
constituées d'un même matériau mais de porosité variant de façon graduelle suivant leur 
épaisseur ou de membranes composées de plusieurs couches de différents matériaux (on parle 
alors de membranes composites). Toutes les membranes de nanofiltration disponibles sur le 
marché possèdent une structure asymétrique car celle-ci permet d’allier une bonne résistance 
mécanique et un débit volumique de perméat élevé. 
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1 Généralités sur les membranes de filtration 




a 


b 


Figure 2 : Représentation schématique des différentes structures membranaires 
(a) Membrane symétrique (b) Membrane asymétrique [Baker, 2004]. 


Une membrane asymétrique est en général constituée d’une couche support 
macroporeuse (assurant une bonne résistance mécanique) et d’une couche active (régissant 
les propriétés de séparation) entre lesquelles s’insèrent une ou plusieurs couches 
intermédiaires dont le rôle est de permettre le dépôt de la couche active et sa tenue en 
pression (figure 3). 


Support macroporeux 


Couchets) intermédiairets'l 


Couche active 


> 


*• 

> 


Fonctions mécaniques 

Fonctions de transition 
Fonctions séparatives 


Figure 3 : Représentation schématique d’une membrane multicouche. 


1.1.3 Géométrie 


Les membranes peuvent se présenter sous forme de feuilles (membranes organiques), 
de plaques (membranes céramiques), de tubes (membranes organiques et céramiques) ou de 
fibres creuses (membranes organiques et céramiques). Selon leur géométrie, elles sont mises 
en œuvre dans les modules suivants : 


16 


1.1 Les différents types de membranes de filtration 


• les modules plans de type filtre-presse (figure 4) ; Ce sont les plus anciens et les plus 
simples : les membranes sont empilées en mille-feuilles séparés par des cadres intermédiaires 
qui assurent la circulation des fluides. 



Figure 4 : Module plan. Schéma d’un montage de type filtre-presse 
[Les cahiers du CFM n°2]. 

• les modules bobinés en spirale (figure 5): une membrane plane est enroulée sur elle- 
même autour d'un tube poreux qui recueille le filtrat. On obtient ainsi un cylindre multi- 
couches où le perméat s'écoule selon un chemin spiralé vers le tube poreux tandis que 
l'alimentation circule axialement dans les canaux. 


Pièce anti-télescopage 
Tube Perforé - Collecte Perméat 


Alimentation 

k 



Bouchon fi ux Rclcntat 


^Support de collecte perméat 
Membrane 
^ Espaceur 
Flux perméat 



Figure 5 : Module bobiné en spirale : schéma d’un module [Les cahiers du CFM n°2]. 
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• les modules tubulaires (figure 6); Ils sont constitués par des faisceaux de tubes de 4 à 
25 mm de diamètre intérieur. Ils sont basés sur une technologie simple, facile d'utilisation et 
de nettoyage, mais ils sont grands consommateurs d'énergie pour une très faible surface 
d'échange par unité de volume (compacité réduite). 



Figure 6 : Module tubulaire [Les cahiers du CFM n°2]. 

• les modules fibres creuses (figure 7); Ils rassemblent un ensemble de fibres creuses de 
diamètre externe inférieur à environ 1 millimètre, en un faisceau. Cette configuration leur 
procure la plus forte densité d'écoulement par module. 



Figure 7 : Module à fibres creuses. Module MICROZA, PALL [Les cahiers du CFM n°2]. 

Les principaux avantages et inconvénients des modules adaptés aux quatre géométries 
de membranes énoncées ci-dessus sont regroupés dans le tableau 3. 
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1.1 Les différents types de membranes de filtration 


Géométrie 

Avantages 

Inconvénients 

Plane 

Possibilité de récupérer le perméat produit 
par chaque élément membranaire 

Nettoyage in situ assez difficile 

Spirale 

Compacité 

Faible coût d’investissement 

Pré-traitement nécessaire 
Nettoyage in situ très difficile 

Tubulaire 

Très peu d’exigences de pré-traitement 
Nettoyage in situ facilité 

Consommation énergétique 
élevée 

Compacité limitée 

Fibre creuse 

Compacité élevée 
Faible consommation énergétique 
Faible coût d’investissement 

Pré-traitement nécessaire en mode 
interne/externe 


Tableau 3 : Avantages et inconvénients des différents modules membranaires 

[Fievet et al., 2009]. 


Plusieurs éléments sont à prendre en compte pour sélectionner une géométrie : 

• le type d’opération de filtration retenue. 

• la compacité recherchée (que l’on peut définir comme le rapport de la surface filtrante 
sur le volume d’occupation au sol de l’installation). 

• la présence de particules qui peut nécessiter un prétraitement dans les cas des modules 
plans ou spirales. 

• la viscosité du fluide : dans le cas des membranes tubulaires ou des fibres creuses, 
plus elle est élevée, plus le diamètre des canaux devra être grand. Une utilisation des fibres 
creuses en mode externe / interne (concentrât à l’extérieur des fibres et perméat à l’intérieur) 
peut également être envisagée. 

Quel que soit le type de module utilisé, celui-ci est composé d’un carter rigide dans 
lequel vient s’insérer la ou les membrane(s), d’une entrée et de une ou deux sorties selon le 
mode de filtration (figure 8). Lorsque la solution arrive perpendiculairement à la surface de la 
membrane, on parle de filtration frontale ; on distingue alors deux courants, la solution 
d’alimentation et le filtrat ou perméat. Dans le second mode de filtration, dite tangentielle, le 
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flux d’alimentation est parallèle à la surface de la membrane. Dans ce cas, on distingue, trois 
courants, la solution d’alimentation, le rétentat ou concentrât et le perméat. Il est à noter que 
la plupart des procédés de filtration mis en œuvre actuellement dans l’industrie utilisent le 
mode tangentiel. 



b) 


Solution 

d'alimentation 



rétentat 


Figure 8 : Principe de la filtration membranaire en a) mode frontal et b) en mode tangentiel 

[Szymczyk et al., 2007]. 


1.2 Propriétés des membranes 
vis-à-vis d’un fluide 


Une membrane peut être caractérisée vis-à-vis d’un fluide à traiter par sa perméabilité 
hydraulique, son seuil de coupure, son taux de rejet, sa sélectivité et sa charge de surface. 


1.2.1 Flux volumique et perméabilité 
hydraulique 

Le flux volumique de perméat est une mesure de volume de perméat par unité de 
temps et par unité de surface membranaire (le flux volumique est donc homogène en 
dimensions à une vitesse). On l’exprime généralement en litre par heure et par mètre carré de 
membrane. Le flux à l’eau pure (J w ) est une caractéristique de la membrane qui est 
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généralement indiquée par le fournisseur. Il obéit à la loi de Darcy qui montre notamment 
que le flux à l’eau est proportionnel à la pression transmembranaire (AP) : 


J 


W 


AP 

TlR m 


eq 1 


où 7] représente la viscosité dynamique de l’eau à la température considérée et R m la 
résistance hydraulique (intrinsèque) de la membrane. 

La perméabilité hydraulique {L P ) représente le flux volumique de perméat par unité de 
pression. On l’exprime en litre par heure, par mètre carré de membrane et par unité de 
pression de travail. C’est une propriété intrinsèque de la membrane qui permet de comparer 
les matériaux entre eux (il est nécessaire de préciser la température à laquelle sont réalisées 
les mesures en raison de l’influence de celle-ci sur la viscosité du fluide). En l’absence de 
soluté (solvant pur) : 


J W =L P AP 


eq 2 


Si l’on considère que la membrane est un milieu poreux, la perméabilité L p peut être 
reliée à la géométrie de pore considérée à partir de l’équation de Hagen-Poiseuille : 


J w = L p ÀP 


8r|(Ax/A k ) 


AP (géométrie cylindrique) 


eq 3 


J w = L p AP = 


3rj(Ax/ A k ) 


AP (géométrie « en feuillets ») 


eq 4 


où r p représente le rayon de pore moyen (pores cylindriques) ou la demi-hauteur moyenne 
(pores « en feuillets ») des pores tandis que Ax/A k représente le rapport de l’épaisseur sur la 
porosité de la membrane. 
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1.2.2 Seuil de coupure 

Les membranes d’ultra- et de nanofiltration sont caractérisées par leur seuil de 
coupure. Il s’agit de la masse molaire de la plus petite molécule arrêtée à 90% (ou à 95% 
selon le fournisseur) par la membrane dans des conditions opératoires données. Il s’exprime 
généralement en Dalton (unité équivalente à des grammes par mole). Le seuil de coupure 
renseigne l’utilisateur sur la plage d’utilisation de la membrane mais sa valeur exacte ne doit 
être considérée qu’à titre indicatif. En effet, les molécules ont des géométries variées et il est 
possible que des espèces de masse molaire inférieure au seuil de coupure soient retenues par 
la membrane, de même que des espèces de masse moléculaire légèrement supérieure la 
traversent. Il est donc important de retenir que le seuil de coupure est défini par rapport à un 
soluté donné. Par ailleurs, l’hétéroporosité (loi de distribution des pores suivant leur 
diamètre) des membranes peut être plus ou moins importante. 

On peut ainsi caractériser une membrane d’ultra- ou de nanofiltration, pour une 
séparation donnée, dans des conditions fixées, par une courbe de variation du taux de rejet 
d’un soluté en fonction de sa masse moléculaire (figure 9). 



Figure 9 : Illustration du principe de détermination du seuil de coupure d’une membrane 
( défini ici pour un taux de rejet de 0,90) [Les Cahiers du CFM n°4]. 
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Cette grandeur n’est pas utilisée dans le cas des membranes de microfiltration pour 
lesquelles les pores sont caractérisés par leur diamètre nominal moyen. Celui-ci représente le 
diamètre de la plus grosse particule sphérique indéformable capable de traverser la 
membrane. 


1.2.3 Taux de rejet 


Le taux de rejet d’une espèce (noté R) est une grandeur sans dimension définie 
comme le pourcentage d’espèces retenues par la membrane. Dans le cas de mélanges 
complexes, un taux de rejet individuel est défini pour chaque type de soluté i. L’accumulation 
éventuelle d’espèces au voisinage de la membrane, ou polarisation de concentration 1 , 
conduit à définir le taux de rejet de deux manières différentes. On distinguera, en effet, le 
taux de rejet observé ( R ii0 bs ) donné par l’équation (5), défini à partir de la concentration de la 
solution d’alimentation (cj~ 8 )et de la concentration du perméat (cf x ) et le taux de rejet 
intrinsèque donné par l’équation (6), où la concentration de la solution d’alimentation a 
été remplacé par la concentration réelle du soluté au voisinage de la surface de la membrane 

(c°). 


R 


i,obs 



eq 5 


Ax 


R ,,int 


= L 


eq 6 


Le taux de rejet observé est une grandeur directement accessible par l’expérience mais 
qui présente l’inconvénient d’être fonction des conditions hydrodynamiques de 
fonctionnement telles que la vitesse de circulation et la géométrie de la lame liquide en 
contact avec la membrane (celle-ci dépend à la fois de la géométrie de la membrane et du 
module membranaire). Le taux de rejet intrinsèque est quant à lui indépendant des conditions 
hydrodynamiques de filtration. Il reflète les performances réelles de séparation d’une 

1 phénomène détaillé par la suite 
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membrane vis-à-vis d’un soluté. C’est une grandeur accessible à partir de la mesure du taux 
de rejet observé si les conditions hydrodynamiques sont connues. 


Polarisation de concentration 


Lorsqu’un gradient de pression transmembranaire est appliqué, solutés et solvant sont 
entraînés par convection vers la membrane. Les solutés, partiellement ou complètement 
retenus, s’accumulent à la surface de celle-ci créant ainsi une surconcentration à l’interface 
membrane / solution externe. Le gradient de concentration engendré entraîne un flux de 
diffusion du soluté, antagoniste, tendant à équilibrer les concentrations. La concentration du 
soluté à l’interface membrane / solution externe augmente jusqu’à ce qu’un état stationnaire 
soit atteint. Le flux de soluté traversant la membrane est alors égal au flux de convection 
dirigé vers la membrane sous l’action du gradient de pression diminué du flux de diffusion se 
déplaçant en sens inverse sous l’action du gradient de concentration. 

La théorie du film [Brites et al., 1991] est le modèle le plus fréquemment utilisé pour 
décrire le phénomène de polarisation de concentration. Le transfert en solution est supposé 
localisé au voisinage de la membrane au sein d’une couche limite d’épaisseur ô (couche de 
polarisation). L’épaisseur dépend essentiellement des conditions hydrodynamiques 
appliquées dans le module et notamment de la vitesse de circulation du fluide à la surface de 
la membrane. 

A partir d'un bilan de matière sur l'épaisseur de la couche de polarisation, le flux de 
soluté transféré à travers la membrane ( J s = J v cf x ) peut s'écrire de la façon suivante (pour un 
soluté non chargé) [Blatt et al., 1970]: 


J s = JvCj (x)-Dj 


dC , O) _ T r Ax 
J J V^i 

dx 


eq 7 


où Jv est le flux volumique du perméat, d(x) est la concentration du soluté en un point x situé 
dans la couche limite et D, son coefficient de diffusion. 


L’intégration de l’équation précédente sur l’épaisseur de la couche limite en tenant 
compte des conditions aux limites c, = c ] 8 pour x = -ô et c i = c® pour x = 0 permet d’écrire: 
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avec 


J v =kln 


f n „Ax y 

Cj ~ c , 

-5 Ax 

V e , -Ci y 


eq 8 


où k désigne le coefficient de transfert de masse. Celui-ci dépend des conditions 
hydrodynamiques du système et des propriétés physico-chimiques de la solution. 

A partir des relations (5), (6) et (8), il est possible d’établir la relation suivante entre le 
taux de rejet observé ( R iM bs ) et le taux de rejet intrinsèque (R Unt ) : 


R Rj.obs ex P (Pe) 

l-R,,*(l-exp(Pe)) 

où Pe représente le nombre de Péclet dans la couche limite et est défini par : 


eq 10 



J v 5 
Di 


eq 11 


Le nombre de Peclet est un nombre adimensionnel qui exprime le rapport entre le flux 
de convection et le flux de diffusion à l’intérieur de la couche de polarisation. 


> Si Pe » 1 : la convection est prépondérante dans la couche limite, la polarisation de 
concentration est favorisée et la rétention diminue. 


> Si Pe « 1 : la diffusion diminue considérablement l’augmentation de concentration à 
la membrane. La polarisation de concentration est négligeable. 


Dans les équations (10) et (11), R it0 b s et J v sont des grandeurs mesurées 
expérimentalement tandis que k doit être estimé à partir de la relation suivante : 
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k _ Dj-Sh 
d h 


eq 12 


où Sh et dh désignent respectivement le nombre de Sherwood et le diamètre hydraulique de la 
section d’écoulement lequel est donné par : 


d h 


4S 

P 


avec S l’aire de la section d’écoulement et P le périmètre mouillé. 


eq 13 


On notera que l’équation (11) reste valable dans le cas d’un électrolyte binaire en 
remplaçant le coefficient de diffusion du soluté neutre D, par le coefficient de diffusion 
global de l’électrolyte D se i qui est défini comme suit : 


Dsel 


(l z l 

+ 

Z 2 1 

) D 1 D 2 

M 

Di+| 

z 2 

D 2 


eq 14 


où D; et D 2 sont respectivement les coefficients de diffusion du cation et de l’anion, et zi et 22 
leur nombre de charge respectif. 


Selon le régime d’écoulement, il existe différentes expressions pour le nombre de 
Sherwood. Celles-ci sont du type [Sherwood, 1934 ; Schweitzer, 1999 ; Bouchoux, 2004]: 

> pour un écoulement laminaire : 


Sh = a Re b 


Sc c 




eq 15 


> pour un écoulement turbulent : 


Sh = a Re b Sc c 


eq 16 


Dans les relations (15) et (16), a , b, c et d sont des coefficients semi-empiriques, L 
représente la longueur de la membrane, Re est le nombre de Reynolds caractérisant le régime 
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d’écoulement et Sc le nombre de Schmidt. Les nombres de Reynolds et de Schmidt sont des 
nombres adimensionnels définis respectivement par : 


Re = 


Sc = 


_ vd h P 


n 

n 


eq 17 


pD ec l 18 

où v représente la vitesse de circulation de la solution d’alimentation et p sa masse 
volumique. 


Diverses corrélations entre les nombres de Sherwood et de Reynolds ont été établies 
dans la littérature selon le type d’écoulement, la géométrie des modules et le type de soluté. 
Quelques exemples sont rassemblés dans le tableau 4. 


Régime d’écoulement 

a 

b 

c 

d 

Conditions 

particulières 

Auteur(s) 

Laminaire 

0,664 

0,5 

0,33 

0,33 

Géométrie plane 
L < 0,0029dhRe 

Grober 

Laminaire 

1,86 

0,33 

0,33 

0,33 

Géométrie plane 
L > 0,0029d h Re 

Graetz - 
Leveque 

Turbulent 

0,04 

0,75 

0,33 

- 

- 

Dittus - Boelter 

Turbulent 

0,023 

0,875 

0,25 

- 

1 <Sc< 1000 

Deissler 

Turbulent 

0,0096 

0,91 

0,35 

- 

Sc > 1000 

Harriot - 
Hamilton 


Tableau 4 : Facteurs semi-empiriques intervenant clans les corrélations entre les nombres 
adimensionnels de Sherwood, de Reynolds et de Schmidt [Bouchoux, 2004]. 
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1.2.4 Sélectivité 


Il s’agit d’une grandeur sans dimension qui traduit l’aptitude d’une membrane à 
rejeter un soluté de type 1 plutôt qu’un soluté de type 2. On définit la sélectivité d’une 
membrane pour un composé 1 par rapport à un composé 2 comme : 


S 


1/2 



eq 19 


où T t désigne le taux de transmission de l’espèce i. 

Comme pour le taux de rejet, il est possible de raisonner sur les grandeurs observées 
(dépendant des conditions hydrodynamiques de filtration) ou sur les grandeurs intrinsèques 
(propres à la membrane). 


1.2.5 Charge de surface 

La plupart des matériaux membranaires, organiques ou inorganiques, acquièrent une 
charge de surface lorsqu’ils sont immergés dans un milieu polaire tel qu’une solution 
aqueuse. Ainsi, lorsque les dimensions des interstices entre les constituants (particules 
d’oxyde ou chaînes de polymère) de la matrice membranaire deviennent comparables à la 
portée des interactions électrostatiques, on conçoit aisément que la sélectivité de la membrane 
vis-à-vis d’un mélange de solutés chargés (ions, polyélectrolytes) ne va pas dépendre 
seulement d’un simple "effet tamis" (c’est-à-dire d’une discrimination en fonction de la taille 
de l’espèce à retenir), mais aussi des effets de charges électriques découlant de la nature 
chimique du matériau membranaire. 

L’origine de la charge de surface d’une membrane peut provenir du caractère acide ou 
basique des groupements de surface (membranes organiques possédant par exemple des 
groupes amide, carboxyle ou sulfone (figure 10)) 
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Polyetlier sulfone sulfoné 



Polyamide aromatique carboxilé réticulé 

Figure 10 : Exemples de polymère entrant dans la composition des membranes organiques 

[Les cahiers du CFM]. 


ou des propriétés amphotères des sites de surface (fonctions hydroxyle à la surface des 
oxydes métalliques constituant les membranes céramiques (figure 11)). L’adsorption 
d’espèces chargées (ions, polyélectrolytes, surfactants ioniques) à la surface du matériau peut 
également être responsable de la charge superficielle. La charge de surface d’une membrane 
dépend donc non seulement de la nature chimique du matériau membranaire mais également 
de son environnement physico-chimique [Les Cahiers du CFM n°4]. 
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Figure 11 : Comportement amphotère d’un oxyde métallique : hydratation et formation de la 
charge électrique en fonction du pH du milieu extérieur [Les cahiers du CFM]. 
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La création de charges à la surface du matériau membranaire entraîne une 
redistribution spatiale des espèces chargées en solution au voisinage de la surface. Les forces 
électrostatiques dues à la charge de surface attirent les espèces de charge opposée et 
repoussent celles de même charge. Le potentiel électrostatique exercé par la charge de surface 
varie progressivement au sein d’une zone appelée « double couche électrique » (DCE) (figure 
12 ). 



Figure 12 : Représentation de la double couche électrique à l’interface solide-solution. 
PHI : Plan d’Helmholtz Interne ; P HE : Plan d’Helmholtz Externe ; CC : Couche compacte ; 
CD : Couche diffuse ; Wo : Potentiel de surface ; Wp : Potentiel au PHI ; % : Potentiel au 
PHE ;C ■' Potentiel zêta [Les Cahiers du CFM n°4]. 
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Le modèle le plus couramment utilisé pour décrire la structure de la DCE ("modèle 
Gouy-Chapman-Stem-Grahame") [Duclert, 1990 ; Chevalier, 1999] suggère que cette 
réorganisation conduit à la formation de deux zones distinctes, la couche compacte (ou 
couche de Stem) et la couche diffuse. 

Couche compacte 

Selon le type d’interaction ion-surface, les ions sont susceptibles de s’approcher plus 
au moins près de la surface. Si les interactions mises en jeu sont de nature purement 
électrostatique, les ions sont qualifiés d’ions indifférents [Hunter, 1981]. Ils conservent leur 
sphère d’hydratation et se positionnent dans un plan fictif parallèle à la surface, appelé Plan 
de Helmholtz Externe (PHE). Ainsi, ce denier correspond à la distance minimale d’approche 
du centre d’un ion solvaté attiré par la surface sous l’effet de la force d’interaction 
électrostatique. 

Dans le cas où les interactions ions / surface ne sont pas uniquement de nature 
électrostatique, les ions sont dits adsorbés spécifiquement [Hunter, 1981]. Ils s’adsorbent 
spontanément même si la surface du matériau n’est pas chargée. Un ion adsorbé 
spécifiquement est susceptible de s’approcher plus près de la surface qu’un ion indifférent en 
perdant une partie ou la totalité de sa sphère d’hydratation. 

Par définition, la couche compacte est une couche adjacente à la surface limitée par le 
PHE. L’épaisseur de la couche compacte est estimée entre 0,5 et 1 nm [Sasidhar et al., 1982 ; 
Hunter, 1993]. 

Couche diffuse 

Le PHE marque la frontière entre la couche compacte et la couche diffuse au sein de 
laquelle les contre-ions (en excès par rapport à la solution de cœur) et les co-ions (déficitaires 
par rapport à la solution de coeur) se répartissent suivant un gradient de concentration jusqu’à 
atteindre les caractéristiques de la solution de coeur. La longueur de Debye (je 1 ) est 
généralement considérée égale à l’épaisseur de la couche diffuse [Lyklema, 1995] (équation 
20). Ce paramètre est très utile pour estimer la portée des interactions électrostatiques. 
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1 Généralités sur les membranes de filtration 


K 


-l 


s , \-l/2 

f 2F 2 I a 


e 0 e ,RT y 


eq 20 


F est la constante de Faraday, 1 la force ionique de la solution, e 0 la permittivité du vide, s r la 
constante diélectrique du solvant, R la constante des gaz parfaits et T la température. 


L’équation (20) montre que l’épaisseur de la couche diffuse diminue avec 
l’augmentation de la force ionique. Ce phénomène est désigné par le terme de compression 
de la DCE [Hunter, 1993]. Le rapport du rayon de pore moyen d’une membrane (r p ) sur la 
longueur de Debye (k 1 ) permet d’estimer la part relative du pore occupée par la couche 
diffuse (comparativement au canal électriquement neutre, c’est-à-dire à la solution de cœur). 
Ce rapport, appelé rayon électrocinétique, indique si les DCE se formant à l’intérieur des 
pores peuvent se développer complètement ou si au contraire elles se chevauchent 
(phénomène de recouvrement de DCE [Hunter, 1981]). 

Ainsi, dans le cas des membranes à petits pores (basse UF et NF), la DCE occupe la 
totalité du pore et la variation radiale du potentiel électrostatique dans le pore devient alors 
négligeable. 


L’interface solide-liquide peut être modélisée par deux condensateurs plans, avec : 


oo = C int OFo-'Pp) 

eq 21 

a P = Cint 0Fp - 'Po) + Cext CPp - 'Ed) 

eq 22 

O d = Cext O^d - 'Rp) 

eq 23 


avec : 


oo + op+Od = 0 


eq 24 


où oo, <Jp et <7,1 représentent respectivement la densité de charge de surface, la densité de 
charge de la couche compacte et la densité de charge de la couche diffuse, Ci n , et C ext les 
capacités intégrales des régions interne et externe de la couche compacte. Wo, ^ et W c \ sont 
respectivement le potentiel de surface, le potentiel au PHI et le potentiel au PHE. 


32 


1.2 Propriétés des membranes vis-à-vis d’un fluide 


Une grandeur généralement utilisée dans les modèles de transport couplant l’équation 
de Nemst-Planck étendue et la théorie de la charge fixe (voir chapitre 2) est la densité de 
charge volumique effective (X) définie comme le nombre de charges fixes par unité de 
volume de pore. Selon la géométrie de pore considérée, elle s’exprime ainsi : 

(pores cylindriques) eq 25 

(pores « en feuillets ») eq 26 

La densité de charge de la couche diffuse (ad) peut être calculée à partir de mesures de 
potentiel zêta en utilisant la théorie de Gouy-Chapman [Hunter, 1981]. 


X = 


2o d 

Frn 


X = 


g d 

Fr n 
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2 Modélisation du transport 
à travers des membranes de 
nanofïltration 


Bien que la nanofiltration soit en plein essor actuellement, il est évident que la 
complexité des mécanismes de transfert à travers des micropores nuit à un développement 
optimal de cette technique séparative sur le plan industriel. Il est donc nécessaire de 
développer des outils suffisamment simples et fiables permettant de comprendre et de prévoir 
les performances de séparation des membranes de nanofiltration. 

Dans ce chapitre, nous rappelons les bases du modèle de transport SEDE (Steric, 
Electric and Dielectric Exclusion) [Szymczyk et al., 2005] puis nous appliquons celui-ci à la 
modélisation du taux de rejet et du potentiel de membrane. 


2.1 Historique des modèles de 
transport 


L’intérêt pratique d’un modèle résulte d’un difficile compromis entre sa simplicité et 
son réalisme. A l’heure actuelle, les propriétés de rétention des membranes de nanofiltration 
n’ont été étudiées qu’à partir d’approches macroscopiques. Parmi les obstacles à l’utilisation 
d’une description microscopique, pourtant plus détaillée et fréquemment utilisée pour l’étude 
du transport dans les canaux ioniques [Mamonov et al., 2003 ; Graf et al., 2004], figurent la 
structure complexe des membranes artificielles et l’échelle de temps caractéristique du 
passage d’un soluté à travers ces matériaux. 

Plusieurs approches (macroscopiques) ont été développées pour décrire le transport de 
matière dans des membranes de filtration. La première est d’ordre phénoménologique et 
découle directement des lois de la thermodynamique des processus irréversibles. Les modèles 
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2 Modélisation du transport à travers des membranes de nanofiltration 


résultants considèrent la membrane comme une boite noire séparant deux compartiments. 
Leur intérêt majeur réside dans le fait qu’aucune information concernant la structure de la 
membrane n’est nécessaire, ce qui leur confère un caractère très généraliste. En s’appuyant 
sur les travaux d’Onsager, Staverman montra que l’ensemble des phénomènes de transport à 
travers une membrane pouvait être décrit, dans un système à n composants, par n(n+ 1)/2 
coefficients phénoménologiques [Onsager, 1931, Staverman, 1952]. Pour un système 
constitué d’un solvant et d’un soluté neutre, Kedem et Katchalsky établirent les relations 
phénoménologiques du flux volumique et du flux molaire de soluté en fonction de trois 
coefficients de transport : le coefficient de réflexion de la membrane, sa perméabilité au 
soluté et sa perméabilité hydraulique [Kedem et al., 1958]. Ces auteurs étendirent ensuite leur 
approche au cas d’une solution contenant un électrolyte binaire [Kedem et al., 1963]. 
L’établissement des équations de flux reposant sur la théorie linéaire de la thermodynamique 
des processus irréversibles, celles-ci ne sont valables que pour de faibles valeurs du flux 
volumique et de la différence de concentration transmembranaire. Cette dernière condition 
n’étant en général pas vérifiée en nanofiltration, Spiegler et Kedem développèrent une 
approche consistant à découper fictivement la membrane en une succession de tranches 
élémentaires, chaque tranche séparant deux éléments fictifs de solution (de concentrations 
virtuelles peu différentes l’une de l’autre) en équilibre thermodynamique avec les faces de la 
tranche considérée [Spiegler et al., 1966]. Les équations établies dans le cadre de la théorie 
linéaire de la thermodynamique des processus irréversibles restent donc applicables mais 
prennent alors un caractère local de même que les coefficients de transport. En intégrant ces 
équations locales sur l’épaisseur de la membrane, Spiegler et Kedem ont établi une équation 
reliant le taux de rejet du soluté au flux volumique, au coefficient de réflexion de la 
membrane et à sa perméabilité au soluté. Le modèle de Spiegler- Kedem a été appliqué avec 
succès dans de nombreux cas et est encore très fréquemment utilisé aujourd’hui. Sa 
principale faiblesse réside dans la procédure d’intégration des équations locales de flux, qui 
suppose que les coefficients de transport locaux sont indépendants des concentrations 
virtuelles. Si cette hypothèse simplificatrice est acceptable dans le cas de membranes neutres, 
elle est plus critiquable dans le cas de membranes chargées séparant des solutés partiellement 

9 

ou totalement ionisés . La notion de charge membranaire étant absente du modèle de 

2 Sauf dans les cas limites de très faibles et de très grand flux volumiques [Lefebvre, 2003]. Ceci s’explique par 
le fait que dans ces deux cas les concentrations locales à l’intérieur de la membrane varient peu avec la position 
axiale dans la membrane. A faible flux, cela est dû à la faible capacité de rétention de la membrane tandis qu’à 
flux élevé, cela provient de l’allure du profil de concentration qui reste “plat” sur l’essentiel de l’épaisseur de la 
membrane et ne chute très brutalement que dans une région proche de la sortie de la membrane. 
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2.1 Historique des modèles de transport 


Spiegler-Kedem, elle fut ensuite introduite dans les expressions du coefficient de réflexion et 
de la perméabilité au soluté en utilisant la théorie de la charge fixe [Hoffer et al., 1967, 
Jitsura et al., 1983]. Nakao et Kimura introduisirent des coefficients hydrodynamiques dans 
les expressions du coefficient de réflexion et de la perméabilité au soluté (modèle SHP 3 4 
[Nakao et al., 1982]) afin de déterminer les paramètres structuraux (rayon de pore et rapport 
épaisseur sur porosité) de membranes d’ultrafiltration à partir du taux de rejet de solutés 
neutres. 


Une autre approche, plus mécanistique, repose sur l’utilisation de l’équation de 
Nemst-Planck étendue (c’est-à-dire incluant le transport du soluté par convection). Cette 
dernière fut proposée par Schlôgl à partir d’une simplification des équations de la 
thermodynamique des processus irréversibles [Schlôgl, 1964]. L’équation de Nernst-Planck 
étendue est utilisée dans le modèle de charge d’espace dont les bases ont été établies par 
Osterle et ses collaborateurs [Morrison et al., 1965 ; Gross et al., 1968 ; Fair et al., 1971]. Ce 
modèle repose également sur l’équation de Navier-Stokes pour la description du flux 
volumique et sur l’équation de Poisson-Boltzmann pour le calcul du profil radial du potentiel 
électrostatique à l’intérieur des pores. En dépit de sa rigueur sur le plan fondamental, le 
modèle de charge d’espace a essentiellement été employé pour étudier les phénomènes 
électrocinétiques à l’intérieur de pores cylindriques et n’a été que très peu appliqué à l’étude 
la rétention d’électrolytes en raison de la complexité des calculs requis [Neogi et al., 1981; 
Wang et al., 1995]. Appliqué à des pores chargés et sélectifs, le modèle de charge d’espace 
est en effet un modèle bidimensionnel tenant compte des variations axiales et radiales de la 
pression, du potentiel électrostatique et des concentrations ioniques à l’intérieur de pores 
cylindriques. Diverses simplifications du modèle de charge d’espace ont été proposées afin 
de diminuer la difficulté et le temps des calculs. Jacazio et al. proposèrent par exemple de 
remplacer la loi de Navier-Stokes par l’équation de Poiseuille pour la description du flux 
volumique [Jacazio et al., 1972] tandis que d’autres auteurs utilisèrent des solutions 
approchées de l’équation de Poisson-Boltzmann [Hijnen et al., 1985 ; Yaroshchuk et al., 
1993], 


3 La théorie de la charge fixe (ou modèle TMS) a été développée par Theorell [Theorell, 1935] et par Meyer et 
Sievers [Meyer et al., 1936]. 

4 Steric Hindrance Pore model. 
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2 Modélisation du transport à travers des membranes de nanofiltration 


Une autre approche permettant de simplifier considérablement les calculs et de les 
étendre aisément au transport de mélanges d’électrolytes découlent des travaux de Schlôgl et 
de Dresner [Schlôgl, 1966 ; Dresner, 1972] 5 . Dans l’approche développée par ces auteurs, 
l’équation étendue de Nernst-Planck est combinée à la théorie de la charge fixe : l’exclusion 
des solutés est régie par des équilibres thermodynamiques aux interfaces (équilibres de 
Donnan) et le transport à l’intérieur de la membrane est décrit par des équations de Nernst- 
Planck étendues couplées à une condition d’électroneutralité locale. La notion de rayon de 
pore est absente et la description du transport est monodimensionnelle. Wang et al. utilisèrent 
cette approche en introduisant les coefficients hydrodynamiques du modèle SHP [Nakao et 
al., 1982] dans l’équation de Nemst-Planck étendue (modèle ES 6 7 8 [Wang et al., 1995]). Dans 

' 7 

la même période, Bowen et Mukhtar proposèrent un modèle équivalent (modèle HM 
[Bowen et al., 1996]). Dans les travaux de ces auteurs, la notion de rayon de pore n’apparaît 
qu’au travers des coefficients hydrodynamiques introduits dans l’équation de Nemst-Planck 
étendue. Bowen et al. ont ensuite introduit la taille des pores dans les équations de partage 

o 

aux interfaces sous forme d’un coefficient de partage stérique (DSPM [Bowen et al., 1997]). 
Appliquer l’approche de Schlôgl et de Dresner à des structures poreuses revient à négliger les 
variations radiales du potentiel électrostatique dans les pores. On parle alors d’approximation 
homogène. Les modèles homogènes sont donc des formes simplifiées du modèle charge 
d’espace dans lesquels l’équation de Poisson-Boltzmann est remplacée par une condition 
d’électroneutralité locale. Palmeri et al. ont défini le domaine de validité de l’approximation 
homogène en terme de rayon de pore et de densité de charge membranaire et montré qu’elle 
pouvait être appliquée raisonnablement à la plupart des membranes de nanofiltration [Palmeri 
et al., 1999]. Tout comme le modèle homogène DSPM proposé par Bowen et al., le modèle 
homogène HET 9 développé par Palmeri et al. considère un mécanisme d’exclusion des 
solutés chargés incluant des effets stériques et électriques [Palmeri et al., 2002]. Parmi les 
différences entre ces deux modèles figure le choix de la grandeur représentative de la taille 
des ions : le rayon cristallin est utilisé dans le modèle HET tandis que le DSPM utilise le 
rayon de Stokes. 


5 Schlôgl étudia le transport à travers des membranes chargées pour des systèmes constitués d’un solvant et d’un 
électrolyte binaire. Dresner étendit l’approche de Schlôgl en traitant le cas de systèmes multicomposants et 
démontra notamment la possibilité d’une rétention totale de certains contre-ions dans les systèmes contenant 
plusieurs types de contre-ions. 

6 Electrostatic and Steric-hindrance model. 

7 Hybrid Model. 

8 Donnan-Steric partitioning Pore Model. 

9 Hindered Electro Transport model. 
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2.1 Historique des modèles de transport 


Yaroshchuk a étudié le phénomène d’exclusion diélectrique généré par la différence 
des constantes diélectriques de la solution et de la membrane (“production” de charges 
images) et a montré qu’il pouvait jouer un rôle primordial en nanofiltration [Yaroshchuk, 
(2000 et 2001)]. Vezzani et Bandini ont utilisé les équations de Yaroshchuk décrivant 
l’interaction diélectrique liée à la “production” de forces images dans des pores en feuillet 
(slit pores) et les ont introduites dans le DSPM [Vezzani et al., 2002], Bowen et Welfoot ont 
également modifié le DSPM en tenant compte de l’exclusion diélectrique mais en considérant 
que celle-ci est régie principalement par la différence d’énergie de Bom d’un ion entre 
l’extérieur et l’intérieur des pores [Bowen et al., 2002]. Szymczyk et Fievet ont récemment 
proposé le modèle homogène SEDE 10 incluant l’effet des forces images (sur la base des 
équations établies par Yaroshchuk pour des pores cylindriques ou en feuillets) et l’effet 
diélectrique de Born (à l’aide d’un modèle de Born modifié) [Szymczyk et al., 2005]. 


10 Steric, Electric and Dielectric Exclusion 
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2 Modélisation du transport à travers des membranes de nanofiltration 


2.2 Modèle homogène SEDE 


Les couches support des membranes n’étant pas sélectives vis-à-vis des ions et des 
molécules organiques de faible masse moléculaire (solutés “cibles” en nanofiltration), elles 
ne sont pas prises en compte dans la description du transport donnée ci-après. La couche 
active de la membrane est décrite par un ensemble de pores droits et parallèles. Nous 
considérerons deux géométries de pore différentes : cylindrique et « en feuillets ». Chaque 
pore est caractérisé par sa longueur Av (représentant l’épaisseur effective, incluant la 
tortuosité, de la couche active) et une grandeur r p représentant le rayon de pore moyen pour 
des pores cylindriques ou la demi-largeur pour des pores « en feuillets » (figure 13) 11 . Afin de 
s’affranchir des effets de bord, nous supposerons Av » r p . 


A.r 


Solution 

concentrée 



r P 

Solution 

| diluée 

i j 

o i 


LAï ^ 


Figure 13 : Représentation schématique d’un pore de la couche active. 


2.2.1 Equations de partage 


D’une manière générale, le potentiel électrochimique d’un ion i ( p i ) de valence z, 
peut être écrit sous la forme suivante : 


11 L’effet de la distribution de la taille des pores [Aimar et al., 1990] n’est donc pas considéré dans le modèle 
proposé. 
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2.2 Modèle homogène SEDE 


Pi = M-f +k B Tln 


'O 

vc°; 


+ z.e'E + W. 


eq 27 


où jil l représente le potentiel chimique standard de l’ion i, ks la constante de Boltzmann, T la 

température, c,- sa concentration, c° la concentration standard, e la charge élémentaire, y le 
potentiel électrostatique et W, l’énergie libre d’interaction incluant toutes les interactions de 
l’ion avec le milieu (énergie de solvatation). 


Dans l’hypothèse d’un équilibre thermodynamique local aux deux interfaces 
membrane / solution, on peut écrire pour un ion i : 



eq 28 


où JJ. , bar représente le potentiel électrochimique à l’intérieur de la membrane. En 
introduisant l’équation (27) dans l’équation (28), on obtient les expressions du coefficient de 
partage de l’espèce i aux interfaces : 


= ex P| 

Ci 


f 


v 


zyA^i 1 ; 11 + A W ; 

k^f 


int "N 


eq 29 


AW; nt =-k B Tln^| nt +k B Tln 


( -int A 

Yi_ 
v in 

\ v ' J 


+ AW +AW 


eq 30 


où le symbole A traduit une variation par rapport à la solution externe, int peut être remplacé 
par 0 ou Ax selon l’interface considérée et A l l((" est appelé potentiel de Donnan. Dans 
l’équation (30) les différents termes à droite du signe égal traduisent respectivement la 
contribution des effets stériques ((f)" 1 ' représente le coefficient de partage stérique), l’effet de 
la variation du coefficient d’activité ionique (y‘ nt ) entre l’intérieur et l’extérieur des pores, 
l’effet diélectrique de Bom (A\\^g om ) et l’effet diélectrique associé à la “production” de 
charges images (AW™^). 
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2 Modélisation du transport à travers des membranes de nanofiltration 


Les conditions d’électroneutralité dans les solutions externes et dans la solution se 
trouvant à l’intérieur des pores s’expriment comme suit : 


I 


Z-C: 


= 0 



= -x 


eq 31 


où X représente la densité de charge volumique effective de la membrane (définie comme le 
nombre de moles de charges fixes par unité de volume de pore et considérée comme 
indépendante de la position axiale dans les pores). 


2.2. 1.1 Effet stérique 


Le coefficient de partage stérique $ apparaissant dans l’équation (30) est défini 
comme le rapport de la section accessible par le soluté sur la section réelle du pore [Ferry, 
1936]. En considérant un soluté sphérique, ce terme traduit le fait que le centre du soluté est 
exclu d’une région, adjacente à la paroi des pores, d’épaisseur égale au rayon du soluté 
(figure 14). 



Aire accessible au soluté 


00 


(a) 


(b) 


Figure 14 : Représentation de l’aire accessible à un soluté sphérique dans un pore 

cylindrique (a) ou en feuillets (b) 


Selon la géométrie de pore considérée, le coefficient de partage stérique s’écrira 

donc : 
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2.2 Modèle homogène SEDE 


(J )"’ 1 = (l - À 1 ) (pores cylindriques) eq 32 

4 ) | 11 ' = (l -À, ) (pores « en feuillets ») e q 33 

avec, 

, _ r i,s 

^1 _ eq 34 

r r 

où r i: s représente le rayon de Stokes du soluté et r p le rayon des pores (pour des pores 
cylindriques) ou la demie hauteur des pores (pour des pores « en feuillet »). 

Le rayon de Stokes (r,-, s) de chacune des espèces est calculé à partir de la relation de 
Stokes-Einstein : 


k g T 

“1 35 

12 

où Dj est le coefficient de diffusion du soluté i à dilution infinie . Les équations (32) et (33) 
sont définies pour des valeurs de À, comprises entre 0 (soluté ponctuel) et 1 (soluté de même 
taille que les pores). 

Si seuls les effets stériques sont considérés dans l’expression de l’excès d’énergie de 
solvatation, le mécanisme d’exclusion aux interfaces résulte alors d’une combinaison d’effets 
électriques (Donnan, 1911) et stériques. L’introduction de l’équation (30) dans l’équation 
(29) conduit alors aux équations de partage utilisées dans les modèles HET [Palmeri et al., 
2002] et DSPM [Bowen et al., 1997)]. 


2.2. 1.2 Coefficients d 9 activité 


Les coefficients d’activité ionique (y”') apparaissant dans l’équation (30) sont 
calculés à l’aide de la loi étendue de Debye-Hückel : 

12 nous négligerons l’influence de la concentration sur le coefficient de diffusion. 
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logyf 


i + a^i-Vt 7 


i°gÿ! nt 


-B int Zl 1 2 3 Vf" 
l + A^i-VÏ^ 


eq 36 


avec 


A mt = F 


e 0 e b RT 


B lnt = 


e F 


4V21nlO(e 0 e b RT)^ 


A = F 


e 0 e p RT 


B"” = 


e 2 F 


4V21nl0(e 0 e RT)^ 


eq 37 
eq 38 


où r, est le rayon de l’ion i , I la force ionique, £b la constante diélectrique de la solution à 
l’extérieur de la membrane et £ p la constante diélectrique de la solution à l’intérieur des pores. 


2.2. 1.3 Effet diélectrique de Born 


En 1920 Bom a proposé un modèle permettant de calculer l'énergie de nature 
électrostatique intervenant dans la solvatation d'un ion sphérique de valence Zi et de rayon 
ionique r, [Bom, 1920]. Dans cette approche macroscopique, le solvant est assimilé à un 
milieu continu homogène de constante diélectrique s r . Le travail nécessaire pour charger un 
ion dans un milieu de constante diélectrique e r est donné par : 


W ; 


i,Bom 


( z i e ) 2 

8ra 0 £ r ri 


eq 39 


Afin de déterminer le travail de charge d'un ion lorsque celui-ci se déplace du vide vers 
un milieu de constant diélectrique e r , c’est-à-dire l’enthalpie libre de solvatation, Born utilisa 
le cycle thermodynamique suivant (figure 15): 


1 . décharger l’ion sphérique dans le vide (le travail nécessaire est égal à wj) 

2. transférer la sphère non chargée du vide dans un milieu de constante diélectrique s r . Le 
travail nécessaire à ce transfert est supposé nul (w„=0) 

3. recharger la sphère dans un milieu de constante diélectrique s r (le travail nécessaire est 
égal à w c ) 
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2.2 Modèle homogène SEDE 



Figure 15 : Cycle thermodynamique utilisé dans le modèle de Born. 


Ainsi selon le modèle de Born, l’enthalpie libre de solvatation est donnée par l’équation : 


AG sq1 = AW; 


i,Born 


( z i e ) 2 

8718 ^ 


— 1 


eq 40 


Si au lieu de considérer le transfert d'un ion du vide vers un milieu de constante 
diélectrique s r , on considère le transfert d'un ion d’un milieu de constant diélectrique ei vers 
un milieu de constant diélectrique e?, le travail requis est donc égal à : 


AW 


i.Born 


( z i e ) 2 

8^ 0 r, 


1 1 


V^2 


u y 


eq 41 


Des études théoriques suggèrent que le confinement d’un liquide [Senapati et al., 
2001] ainsi que le champ électrique généré par une surface chargée [Blaak et al., 2006] 
entraînent une réduction de la polarisabilité des molécules du solvant et donc une diminution 
de sa constante diélectrique. Appliqué aux systèmes membranaires, le modèle de Born peut 
être utilisé pour calculer le travail nécessaire pour transférer un ion i depuis la solution 
externe de constante diélectrique Ch jusqu’à l’intérieur des pores à l’intérieur desquels la 
constante diélectrique de la solution vaut s p . Cependant des mesures de l’enthalpie de 
solvatation réalisées sur de nombreux solides ioniques montrent que celle-ci est surestimée 
par le modèle de Born, notamment pour les cations [Latimer et al., 1939]. Rashin et Honig 
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ont suggéré de remplacer le rayon de l’ion par le rayon de cavité [Rashin et al., 1985]. Le 
rayon de cavité (r icav ) représente la distance entre le centre de l’ion et le point où la constante 
diélectrique devient différente de celle du vide. Ils ont déterminé cette distance à partir de 
l’analyse des distributions de densité électronique dans des cristaux ioniques. Rashin et 
Honig conclurent que le rayon de cavité est proche du rayon ionique pour les anions et 
proche du rayon covalent pour les cations (un accord optimal avec l’expérience étant obtenu 
en augmentant les rayons précédents de 7 %). Szymczyk et Fievet ont utilisé l’approche de 
Rashin et Honig afin de calculer AW; m g om dans le modèle SEDE [Szymczyk et al., 2005)]. La 

variation de l’énergie de solvatation ( AW™^ ) s’écrit alors : 


AW- 


i,Bom 


( z i e ) 2 

S^oLcav 


1 1 


V P 


eq 42 


Cette équation suppose que i) le rayon de cavité ne dépend ni de la concentration ni de la 
constante diélectrique, ii) la constante diélectrique est identique en tout point du pore. 


2.2. 1.4 Effet diélectrique dû aux charges images 


Au début du XX eme siècle, Heydweiller a montré expérimentalement qu’une solution 
contenant un sel inorganique possède une tension de surface supérieure à celle du solvant pur 
[Heydweiller, 1910]. D’un point de vue thermodynamique, cela traduit une déficience en 
soluté au voisinage de l’interface. Cette zone de déplétion ne pouvant pas être décrite par 
l’approche classique de Poisson-Boltzmann, Wagner proposa la première théorie des forces 
images [Wagner, 1924]. En raison de la discontinuité diélectrique à l’interface liquide / 
vapeur, le champ électrique généré par les ions polarise l’interface et les charges de 
polarisation induites interagissent de manière répulsive avec ces mêmes ions. Le terme 
« charge image » provient du fait que, dans certains cas, l’interaction générée est identique à 
celle qui se produirait entre un ion situé à une certaine distance de l’interface et une charge 
virtuelle (appelée image de l’ion) qui serait située dans la phase vapeur à égale distance de 
l’interface. Dans la théorie de Wagner, le potentiel d’interaction dû à la production de charges 
images à l’interface entre deux milieux de constantes diélectriques s\ et £2 est égal à 
[Wagner, 1924 ; Manciu et al., 2003] : 
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AW, 


Wagner 


(x) 


\ £ 2 + £ l J 


(Z,g ) 2 

16 7T£q£ 2 X 


exp 


2x 

r-'O) 


eq 43 


où x représente la distance par rapport à l’interface et k la longueur d’écrantage (longueur 
de Debye). 

Comme le montre l’équation (43), tous les ions de même valence sont caractérisés par 
le même potentiel d’interaction, ce qui n’est pas forcément vérifié expérimentalement. C’est 
la raison pour laquelle Onsager et Samaras introduisirent la taille de l’ion (r,j dans la théorie 
de Wagner, [Onsager et al., 1934 ; Manciu et al., 2003] : 


A W os (x) = 


r \ ex P 

£ 2 ~ £ \ 


\ £ 2 £ \ J 


2 Tj_ 

\ K ' J 


{z t ef 


1 + 


2 r 


16 7T£ q £ 2 X 


exp 


2x ) 

A - 1 J 


K 


eq 44 


La théorie d’ Onsager et Samaras donne des résultats satisfaisant à faibles 
concentrations [Long et al., 1942 ; Passoth, 1959] mais elle sous-estime la tension de surface 
pour les fortes concentrations (c>0,2 M) [Randles, 1963]. Une des raisons est que la théorie 
développée par Onsager et Samaras suppose que la longueur d’écrantage (K 1 ) est 
indépendante de la distance à l’interface. Bostrôm et al. ont par ailleurs suggéré d’ajouter aux 
forces images, une force de dispersion laquelle devient prépondérante à forte concentration 
(les charges images étant écrantées par l’atmosphère ionique) [Bostrôm et al., 2001]. 


Afin d’appliquer la théorie des charges images aux interfaces membranaires, 
Yaroshchuk modifia l’approche d’Onsager et Samaras en considérant toujours que la 
longueur d’écrantage est constante mais en la calculant sur la base de la concentration 
moyenne dans les pores (et non pas sur la base de la concentration de bulk comme c’est le cas 
dans la théorie d’Onsager et Samaras) [Yaroshchuk, 2000]. Yaroshchuk a établi des 
expressions approchées de l’énergie d’interaction résultant des forces images dans des pores 
chargés de différentes géométries [Yaroshchuk, 2000 et 2001]. Ses équations tiennent compte 
de l’écrantage de l’interaction diélectrique par la charge de surface et par l’atmosphère 
ionique autour d’un ion donné. Yaroshchuk a montré que l’énergie d’interaction peut être 
écrite comme la somme de deux termes dont l’un est fonction de la longueur d’écrantage 
mais pas de la distance par rapport à la surface (ce terme représente l’énergie d’interaction 
pour un ion situé au centre du pore) et dont l’autre dépend de la position radiale mais pas de 
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la longueur d’écrantage. En négligeant le second terme devant le premier, Yaroshchuk 
obtient des expressions approchées de l’excès d’énergie de solvatation dû à l’interaction d’un 
ion avec les charges de polarisation induites à la surface d’un pore chargé. Dans le cadre du 
modèle SEDE celles-ci s’écrivent : 


avec 


AW. ant 

i, image 


2a, r K 0 (k')K 1 (v)-P(k')K 0 (v)K 1 (k') 


r^o 

71 


(v)K 0 (k’) + (k’) 
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Wârt i i ÿ; ^ 


7 eA x E mt A 

^ - — AW' *5 - AW' int 

i,Bom i,im 


v k B T 
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eq 45 


eq 46 


eq 47 

eq 48 
eq 49 

eq 50 


Dans ces équations, k’ représente le vecteur d’onde (espace de Fourier), 7 0 , 1\, K {) et 
K\ les fonctions de Bessel modifiées, Ad^'J" le potentiel de Donnan, AWj’” 1 l’effet 

diélectrique associé à la “production” de forces images normalisé par k B T et AW i Born 
l’excès d’énergie de solvatation dû à l’effet de Born normalisé par k B T . 


Comme le montrent les équations (42) et (45) à (47), l’interaction diélectrique liée à 
l’effet de Bom ou aux forces images, varie avec le carré de la charge d’un ion. Les effets 
diélectriques auront donc tendance à exclure aussi bien les cations que les anions. Cela 
représente une différence majeure avec l’exclusion de Donnan pour laquelle certains ions (co- 
ions) sont exclus de la membrane tandis que les autres (contre-ions) sont attirés à l’intérieur 
des pores. 


48 


2.2 Modèle homogène SEDE 


2.2.2 Equation de transport 

Toutes les grandeurs (flux, vitesses, concentrations, potentiel électrostatique et 
énergies d’interaction) sont définies comme des fonctions moyennées radialement [Cwirko et 
al., 1989], les déviations spatiales par rapport à la moyenne étant négligées (approximation 
homogène). 

L’équation étendue de Nernst-Planck [Schlôgl, 1964] est à la base des modèles de 
transport actuels (ES, DSPM, HET, SEDE). Elle décrit le transport d’un soluté en termes de 
diffusion sous l’action d’une différence de concentration, de migration sous l’action d’un 
champ électrique et de convection liée au flux volumique. L’équation étendue de Nernst- 
Planck a été modifiée à l’aide de coefficients hydrodynamiques Kij et Kj X rendant compte 
respectivement de l’effet des parois des pores sur le transport diffusif et convectif [Anderson 
et al., 1974 ; Deen et al., 1980] et prend alors la forme (en utilisant des unités molaires pour 
le potentiel électrochimique) : 


de z.c.K. .D.E dT 

j, =-K id D i — ^ — — — — + K c V 

' dx RT dx 


où jj est le flux molaire local de l’ion i et V est la vitesse du fluide. 


eq 51 


Les coefficients hydrodynamiques dépendent du rapport de taille entre le soluté et le 
pore (À,/) mais également de la position radiale à l’intérieur du pore. L’effet de la taille finie 
des pores sur la diffusion et la convection peut néanmoins être estimé raisonnablement en ne 
considérant que les valeurs de K i>( i et K ix au centre des pores [Deen, 1987 ; Bowen et al., 1997 
; Bungay et al., 1973]. 


Dans le cas de pores en feuillets, les expressions approchées de K hd et K ix utilisées 
sont les suivantes [Deen, 1987] : 


13 L’introduction de facteurs hydrodynamiques pour décrire la diminution de la mobilité de solutés dans des 
structures poreuses a été proposée par Pappenheimer et al. [Pappenheimer et al., 1951] et Renkin [Renkin 
,1954]. Les effets stériques sont également introduits à travers ces coefficients de gène à la diffusion et à la 
convection 
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K id = 1 -1.004/1,. + 0.4 184 3 + 0.2 U,- 4 - 0.1 69/1, 5 


A', , 


3-^- 


! 2 ) 


, 2 ^ 

>4 

v y 


eq 52 
eq 53 


Dans le cas de pores cylindriques [Bungay et al., 1973] les expressions approchées de 
Km et Ki' C sont les suivantes : 


4/ = 


6 n 
K, 


l,t 


K u = 


( 2 - 4 *-,, 

2K : , 


eq 54 
eq 55 


et où A), et A] v sont définis par : 








n—\ 


n=l 


+ 2 a »+ 3 ^" 
n=0 

71=0 


eq 56 
eq 57 


avec «i = -73/60, «2= 77.293/50.400, <23= -22.5083, «4= -5.6117, a 5 = -0.3363, «6 = -1.216, «7 
= 1.647, b] = 7/60, b 2 = -2.227/50.400, £3= 4.0180, Z> 4 = -3.9788, fc 5 = -1.9215, Z> 6 = 4.392, £7 = 
5.006. (f)i et A; sont donnés par les équations (32) à (34). 


2.2.3 Taux de rejet ionique 

Le processus de séparation est modélisé par la succession d’étapes suivante : (i) tout 
d’abord, les solutés se répartissent à l’interface membrane / rétentat selon les équations de 
partage décrites précédemment ; (ii) le transport des solutés à l’intérieur de la couche active 
est ensuite décrit par l’équation de Nemst-Planck étendue ; (iii) finalement, les solutés se 
répartissent à l’interface membrane / perméat (selon les mêmes lois de distribution qu’à 
l’interface membrane / solution d’alimentation). 
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La chute de pression dans les couches support n’est pas toujours négligeable par rapport à la 
différence de pression appliquée à travers la membrane et est en général difficile à évaluer. 
Afin de s’affranchir de cette difficulté, il est préférable d’utiliser le flux volumique comme 
variable indépendante plutôt que la différence de pression. En effet, à l’état stationnaire, le 
flux volumique est constant et identique dans chacune des couches traversées par le fluide. Le 
flux molaire d’un soluté (défini par unité d’aire de pore) est relié au flux volumique Jy 
(grandeur mesurable définie par unité d’aire de membrane) comme suit : 

Ax 

Ax J . , C- 

j, = c i V = eq 58 

A 

où Ak représente la porosité de la couche active. 

En introduisant l’équation (58) dans l’équation (51) et en supposant que les variations 
de l’excès d’énergie de solvatation sont négligeables à l’intérieur des pores, on obtient 
l’expression du gradient de concentration d’un ion i dans la couche active : 


= h ( K ï,_ c *.)_ï£ELë! 

dx K jd D i „A k " 1 RT dx 


eq 59 


L’équation d’électroneutralité (eq 31) et l’équation (59) permettent d’exprimer le 
gradient de potentiel électrique à travers les pores sous la forme suivante : 


dAF 

dx 


Z 


z. J, 


T 1 K id D ioo A k 


(Ki,c,-c“) 


F y 2 

/ c i z : 

RT ; 


eq 60 


Le taux de rejet intrinsèque d’un ion i (/? rmt ) est défini comme la fraction de cet ion 
retenue par la membrane à l’état stationnaire: 


Ai 


Run t =F 


C : 


eq 61 
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2.3 Application du modèle SEDE 
à l’étude du potentiel de 
membrane 


L’étude du transport d’un soluté à travers une membrane sous l’action d’un gradient 
de concentration permet d’obtenir des informations relatives aux propriétés électriques de 
surface du matériau membranaire [Yamauchi et al., 1989 ; Asaka, 1990 ; Takagi et al., 1996]. 
Lorsqu’une membrane neutre sépare deux compartiments contenant le même électrolyte à la 
même pression hydrostatique mais à des concentrations différentes, le transport des ions à 
l’instant initial résulte de la diffusion des cations et des anions en fonction de leurs gradients 
de concentration. Comme les mobilités des ions sont généralement différentes, il se crée un 
flux de charges d’un compartiment à l’autre. Ce flux induit un champ électrique dans la zone 
de diffusion (établissement d’un potentiel de diffusion ou potentiel de jonction) qui s’oppose 
au déplacement des espèces les plus mobiles et au contraire accélère les plus lentes jusqu'à ce 
que leurs flux soient égaux. Ainsi, après un régime transitoire très bref, il s’établit un régime 
stationnaire de diffusion pendant lequel le flux ionique est électriquement neutre (le régime 
n’est réellement stationnaire que si les concentrations dans les deux compartiments sont 
maintenues constantes dans le temps) [Brun, 1989]. Dans le cas d’une membrane chargée, les 
interactions électrostatiques peuvent influencer de manière non négligeable le transport des 
espèces ioniques qui résulte alors d’un phénomène d’exclusion - diffusion [Westermann- 
Clark et al., 1986]. La différence de potentiel s’établissant de part et d’autre de la membrane 
est alors appelée potentiel de membrane. 

Le transport sous l’action d’un gradient de concentration peut être abordé de 
différentes manières. Le potentiel de membrane a été étudié par une approche 
thermodynamique par Henderson basée sur les travaux de Planck. La relation de Planck- 
Henderson, largement utilisée dans la littérature, permet d’étudier les propriétés 
électrocinétiques des membranes [Szymczyk et al., 1998 ; Benavente et al., 2001]. Szymczyk 
et al. ont étudié l’influence de la concentration et du pH sur le nombre de transport d’ions 
dans une membrane céramique microporeuse séparant des électrolytes binaires symétriques 
et asymétriques [Szymczyk et al., 1998]. D’autres études ont été réalisées au moyen du 
modèle de charge d’espace [Fievet et al., 1999 ; Shang et al., 2006]. Fievet et al. ont étudié 
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théoriquement l’effet de la charge et de la nature du sel sur le potentiel de membrane. Le 
potentiel de membrane a également été étudié à l’aide du modèle homogène dans le cas de 
membranes poreuses [Takagi et al., 1996 et 2001] et de membranes échangeuses d’ions 
[Kimura et al., 2005]. 

Dans la suite de ce travail, le modèle SEDE est appliqué à l’étude du potentiel de 
membrane afin d’étudier l’influence des effets stérique, électrique et diélectriques sur le 
phénomène d’exclusion-diffusion. D’après l’approche TMS [Theorell, 1935 ; Meyer et al. , 

1936], le potentiel de membrane ( A'E,,,) peut être décomposé en deux contributions : 

> la différence des potentiels de Donnan aux interfaces ( A'E ^ and A'F ^ : Av and 0 font 
respectivement référence à l’interface membrane / solution la plus diluée et à l’interface 
membrane / solution la plus concentrée) 

> le potentiel de diffusion à l’intérieur de la membrane ( Ah / clilt ) 


Le potentiel de membrane peut donc s’écrire : 


A*? 

m 


= A’F“ -A^+A^ 


eq 62 


2.3.1 Potentiel de Donnan 


Afin de faire apparaître le potentiel de Donnan, les équations (29-30) peuvent être 
réécrite de la manière suivante : 


A xp int = XJ/ mt _ XJ/ .nt = _ËjlL ln 

Z : e 


f int „ int N 

K , C , 

—int 

V c i j 


eq 63 


avec 


< = 4>“ exp (- 4W*. )exp (- AW», ) 


eq 64 
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2.3.2 Potentiel de diffusion 


En l’absence de flux volumique (on suppose que la différence de concentration est 
suffisamment faible pour que la différence de pression osmotique de part et d’autre de la 
membrane n’engendre pas de flux volumique) l’équation (51) prend la forme 14 : 


Ji 


-K; D. 


f - 

dci 

dx 


+ 


Fz^i d\|/ 
RT dx 


eq 65 


v y 

A l’état stationnaire, aucun courant ne traverse la membrane ce qui se traduit par la relation 
suivante : 


Z Fz J. =0 

i 


eq 66 


En introduisant l’équation (65) dans la condition de courant nul (eq 66), on peut écrire : 


dT_ RT 
dx F 


I K 11 - % 

ZKiD.zfCi 


eq 67 


En considérant l’approximation Henderson le profil de concentration dans le pore peut être 
écrit : 


— Ax —0 

- , . Ci - Ci -0 

Ci (x ) = X + Ci 

Ax 


eq 68 


En introduisant l’équation (68) dans l’équation (68) puis en intégrant celle-ci sur la longueur 
du pore, on obtient l’expression suivante du potentiel de diffusion (A'Ediff) : 
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eq 69 


14 d car le terme convectif est négligé 
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2.3.3 Potentiel de membrane 

En introduisant les équations (63) et (69) dans l’équation (62), on obtient l’expression 
suivante du potentiel de membrane [Matsumoto et al., 1998 ; Takagi et al., 2001] : 
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eq 70 
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3 Recueil de publications 


La complexité des mécanismes de transfert en nanofiltration (NF) nuit à un 
développement optimal de ce procédé sur le plan industriel. C’est pourquoi les recherches 
visant à relier les caractéristiques d’un nano filtre à ses propriétés de transfert et à mettre au 
point des outils fiables de modélisation du transport, s’intensifient depuis quelques d’années. 

Un tel modèle a été récemment développé par Szymczyk et Fievet [Szymczyk et al., 
2005]. Celui-ci se différencie des autres modèles par la prise en compte simultanée de trois 
types d’interaction : l’exclusion stérique, l’exclusion de Donnan et les effets diélectriques 
incluant à la fois les effets de Bom et les forces images. La validation de ce modèle passe 
nécessairement par une comparaison des taux de rejet calculés à ceux obtenus 
expérimentalement. Or, la détermination expérimentale de certains paramètres du modèle, en 
particulier la densité de charge volumique à l’intérieur des pores et la constante diélectrique 
dans les pores, est extrêmement difficile, voire impossible. 

Afin de pouvoir déterminer ces paramètres, il m’a été proposé d’étendre l’application 
du modèle SEDE à l’étude du potentiel de membrane. Ce travail, à la fois théorique et 
expérimental, a donné lieu à 3 publications. 


La première publication est consacrée à l’étude théorique de l’influence des effets 
stérique, électrique et diélectriques sur le potentiel de membrane dans le cas d’un mélange 
binaire (2 ions). Le rapport des concentrations de part et d’autre de la membrane a été fixé à 
2. Les calculs ont été effectués en considérant des membranes neutres puis chargées et un 
électrolyte asymétrique, le chlorure de calcium, à différentes concentrations. Dans le 
formalisme du modèle SEDE, le potentiel de membrane est la somme de deux contributions : 
la différence des potentiels de Donnan aux interfaces membrane / solutions externes et le 
potentiel de diffusion à l’intérieur de la membrane. L’étude révèle que la mesure du potentiel 
de membrane à concentration élevée pourrait constituer une méthode alternative pour 
déterminer la taille moyenne des pores d’une membrane de NF. 
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3 Recueil de publication 


Dans la seconde publication est proposée une analyse simultanée de deux grandeurs 
caractéristiques du transport : le potentiel de membrane et le taux de rejet ionique. L’objectif 
est de caractériser la membrane en termes de densité de charge volumique et de constante 
diélectrique à l’intérieur des pores et ainsi de tester la cohérence globale du modèle SEDE. 
Les expériences ont été menées avec une membrane de 2 eme génération en polyamide et des 
solutions de chlorure de calcium à différentes concentrations (0.30 à 11 mol.m' ). Les 
résultats obtenus démontrent la cohérence globale du modèle SEDE ainsi que l’importance de 
phénomènes d’exclusion autres que stérique et électrique, probablement diélectrique. 


La troisième publication fait suite à la première. Elle concerne l’étude théorique du 
potentiel de membrane s’établissant de part et d’autre d’une membrane de NF séparant deux 
solutions ternaires (3 ions). Le rapport des concentrations de part et d’autre de la membrane 
est fixé à 2 pour chaque ion. Les calculs sont effectués avec une membrane chargée 
négativement et des mélanges constitués de chlorure de potassium et de chlorure de calcium. 
Il ressort de cette étude qu’aux fortes concentrations, le potentiel de membrane tend vers une 
constante dépendant de la composition du mélange, de la taille et de la constante diélectrique 
dans les pores. Puisque la mesure du potentiel de membrane avec des solutions binaires 
concentrées permet de déterminer la taille des pores (publication n°l), les mesures réalisées 
avec des mélanges ternaires permettraient alors d’estimer la constante diélectrique à 
l’intérieur des pores (sans avoir besoin d’effectuer des mesures supplémentaires de taux de 
rejet d’ions). 
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Publication 1 : Influence of Steric, 
Electric, and Dielectric Effects on 
Membrane Potential 


Y. Lanteri, A. Szymczyk et P. Fievet 
Langmuir 2008 , 24, 7955-7962 
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Publication 1 : Influence of Steric, Electric, and dielectric effect on membrane 
potential 
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Publication 2 : Evaluation of the 
steric, electric, and dielectric 
exclusion model on the basis of sait 
rejection rate and membrane 
potential measurements 


Y. Lanteri, P. Fievet et A. Szymczyk 
Journal of Colloid and Interface Science 331 (2009) 148-155 
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Membrane potential in multi-ionic mixtures 


1 2 1 
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Abstract 


The membrane potential arising through charged porous membranes separating 
electrolyte mixtures at identical hydrostatic pressures but different concentrations was 
investigated theoretically. Transport phenomena were described through the Nemst-Planck 
formalism while ion partitioning at the membrane / solution interfaces was accounted for by 
means of modified Donnan équations including steric and dielectric effects. The high- 
concentration limit of the membrane potential dépends on the mixture composition and on the 
pore size as well. A spécifie feature of membrane potential in multi-ionic Systems is its 
dependence at high-concentration with the effective dielectric constant of the solution 
confined inside pores. Indeed, the effective dielectric constant inside pores does not affect the 
high-concentration limit of the membrane potential in the case of single-salt solutions. The 
low-concentration limit of the membrane potential is independent of the mixture 
composition, of the effective dielectric constant inside pores and of the pore radius but it is 
ruled by counterions with the highest charge number. 

The membrane potential measured at high sait concentration with single-salt solutions and 
electrolyte mixtures could be used to détermine the pore size and the effective dielectric 
constant inside pores, respectively. This may constitute an alternative way for membrane 
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characterization with the advantage of avoiding the need of additional rejection rate 
measurements. 
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1. Introduction 


Investigating ion transport induced by a concentration gradient has proved to be an 

1 7 

efficient way to characterize the surface electrical properties of ion-exchange ' and porous 
membranes like microfiltration 8 9 , ultrafiltration 10 " 12 and nanofiltration 13 " 16 membranes. 

Theoretical studies on the membrane potential phenomenon hâve been performed from either 
the Teorell-Meyer-Sievers (TMS) model ’ or the space charge (SC) model . Westermann- 
Clark and Anderson compared membrane potentials predicted by these two approaches in the 
case of 1-1 electrolytes. Their findings showed that the TMS model and the much more 
complex SC model led to similar results for sufficiently low surface charge densities. Shang 
et al. showed that TMS and SC models coincide with each other for slightly charged 
membranes separating symmetric or asymmetric electrolytes provided the pore radius is less 
than 5.0 nm. Otherwise, the TMS model (which neglects radial variations of both electrical 
potential and ion concentrations) was shown to overestimate the membrane potential. 

Steric hindrance is known to play a rôle in the séparation of solutés by nanofiltration 

22 23 24 

membranes and is taken into account in modem transport models. ’ Cervera et al. and 

oc 

Lanteri et al." investigated the membrane potential phenomena with binary electrolytes by 
including ion size effects in the SC model and the TMS model, respectively. They concluded 
that ion size effects lead to an increase in the membrane potential. Lanteri el al. also 
investigated the influence of dielectric effects on the membrane potential arising through a 
membrane separating two solutions of the same single electrolyte. For charged membranes, it 
was shown that both steric and dielectric effects increase the membrane potential at a given 
concentration but the high-concentration limit of the membrane potential is affected only by 
steric effects. 

Ail the above-mentioned works were performed by considering binary electrolytes. In the 

oc 

présent work, we extend our previous analysis - including both steric and dielectric effects in 
the case of the membrane potential arising through charged membranes separating multi- 
ionic mixtures. Spécial attention is paid to the low and high-concentration limits of the 
membrane potential in temary mixture made of a symmetric electrolyte (KC1) and an 
asymmetric one (CaCL) with a common anion. 
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2 Theoretical background 


The electrical potential différence arising through a membrane that séparâtes two 
solutions of the same electrolyte(s) at the same température and hydrostatic pressure but at 
different concentrations is termed membrane potential. 1 The current work deals with the 
theoretical analysis of the membrane potential phenomenon occurring when a charged 
membrane séparâtes multi-ionic mixtures. 

Within the scope of the TMS model (and related models like the SEDE model), the 
membrane potential ( A x ï ) m ) can be split in two components, namely the différence between 

the Donnan potentials at each interface ( AT^ and A*T D : superscripts Ax and 0 refer to the 
interface between the membrane and the most diluted solution and the interface between the 
membrane and the most concentrated solution, respectively) and the diffusion potential 

( A ) arising through the membrane pores 1 


AT( 11 =AT' D Ax -AT'°+A'ï' dlff (D 

We limit ourselves to the situation where diffusion drives ail ions towards the same direction, 
i.e. from 0 towards Dx. 

2.1 Donnan potential 


At each interface the partitioning coefficient of an ion i is defined as the ratio of its 
concentration just inside the membrane pores, cT , to its concentration just outside the 
membrane, c ( int (superscript “int” stands for 0 or Av depending on the interface that is 
considered). Within the scope of the SEDE model the partitioning coefficient accounts for 
steric hindrance, the Donnan exclusion and dielectric exclusion 26 
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— int 

C i int 

— = AT, exp 

int 1 r 


Z, ' 

k B T , 


( 2 ) 


where z.i is the charge number of ion i, e the elementary charge, k B the Boltzmann constant, T 
the absolute température and te!" 1 is the non-Donnan partitioning coefficient defined as 


( 3 ) 


where </>'" represents the steric partitioning coefficient for ion i (which is defined as the ratio 

of the available section for ion i to the pore cross section), A Wjj) l nrn and A W i ^j mge dénoté 

the increase in the interaction energy due to Born dielectric effect and image forces, 
respectively (the prime symbol indicates that both terms are scaled on k B T). 


Considering cylindrical pores, the steric coefficient is given by 



( 4 ) 


with A.= r.l r , r. is the Stokes radius of ion i and r P the pore radius. 


The work of charge transfer from extemal solutions characterized by a dielectric constant ci, 
to a medium with the dielectric constant a p is given by 


84 


Publication 3 : Membrane potential in multi-ionic mixture 


AVE 


int 

i, Boni 


(Zj eY 


f \ 

1 1 


^O k B Tr i,cav 


V £ P 


'b ) 


( 5 ) 


where £ is the vacuum permittivity and r the radius of the cavity formed by the ion i in 

0 i , cav 

27 

the solvent." Within the scope of the SEDE model, e p corresponds to the effective dielectric 
constant of the solution inside pores and £b to the dielectric constant of the extemal bulk 
solution. 


The dielectric exclusion resulting from image forces lias been reviewed recently by 
Yaroshchuk who derived approximate expressions for the interaction energy due to image 
forces. These équations account for the screening of the interaction between the ion and the 
polarized interface b y the membrane fixed charge and the ionic atmosphère as well. Within 
the scope of the SEDE model, these équations read as follows in the case of cylindrical 
pores 26 


AW int 

i ,image 


2 ^ j K'(k)K, (v)-J<Jc)KJy)K, (Je ) ^ 
* • I : (v)K,(k) + p(k)Wv)K : (k) 


(6) 


where / Q , , K Q , K j are the modified Bessel functions, k the wave vector and 


q,-g ) 2 

8 7± B Tr p £ 0 £ p 


( 7 ) 


V = 


\k 2 + 


int 




( 8 ) 
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with 


M 


int __ A Jnt / 


Z- 


z?crrexp ( -^ 

k B T 




21 


int 


(9) 


with 


K 


int 


( r 1/2 

£ 0 £ P RT 

2 F 2 / int 

v y 


( 10 ) 


where F is the Faraday constant, R is the idéal gas constant and 7 int is the ionic strength 
defined as 


^int 



(H) 


P 


int 






J 


( 12 ) 


where £ is the dielectric constant of the membrane material. 

m 


Rearranging eq 2 and considering c° = 2cf* yields the following expression for the 
différence in the Donnan potentials appearing in eq 1 
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A'F 


Ax 


D 


n 

■A'î'g = -^—\n 


Zi e 


^ 1 


.-Ax , 
2 C i JC: 


.0 

i y 


(13) 


2.2 Diffusion potential 


The Nemst-Planck équation (NPE) describes ion transport through membrane pores in terms 
of diffusion and electromigration. The volume flux resulting from the osmotic pressure 
gradient is assumed to be negligible and is disregarded throughout this work. Since pore 
length of synthetic membranes is much larger that their diameter, ion molar fluxes (/,-) can be 
considered one-dimensional and the NPE reads as follows 


j, =~K,D , 


d Ci F z 

h — C i 

dx RT 


d^ ' 
dx 


(14) 


where D is the diffusion coefficient at infinité dilution, c, and 'P are the local ion 

i 

concentration and electrostatic potential inside pores, respectively and K; is the hindrance 
diffusion coefficient accounting for the effect of finite ion and pore sizes. For cylindrical 
pores, the hindrance diffusion factor takes the approximate form 29 


K; = 


6n 


(15) 


K- 


i,t 


-7i 2 V2(l-À.) 
4 1 


-5/2 


1+ Z a n (l-\) n 

n=l 


4 

+ Z 

n=0 


^ t n 
a n+3^i 


(16) 


with a, =- 73/60, a 2 =77.293/50.400, a 3 =-22.5083, « 4 =-5.6117, « 5 =-0.3363, a 6 =- 1.216, 
a 7 1.647. 
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At steady State, the ion fluxes are coupled via the zéro net electrical current flowing through 
the membrane pores 


TFz i j i = 0 


(17) 


The electro-neutrality condition within the membrane pores reads as folio w s 



(18) 


where X is the membrane volume charge density. 


Introducing eq 14 into eq 17 yields 


cl T 
dx 


k B T 

e 


y kd z 

/ J l l l 


- \ 
d Ci 

dx 


'Z K iD l zîc t 


(19) 


T 

Assuming constant concentration gradients we can write 


c> U) 


Ci - Ci 
X + 

Ax 


— 0 

C i 


( 20 ) 


where A x is the pore length. 


Rearranging eq 19 and eq 20 and integrating between pore ends leads to the folio wing 
expression for the diffusion potential for any number of ions 
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diff 


k B T 

( -V-I —Aï -0 V 

Y J K i D,z i (c i -a) 

i 

ln 

f' VI ,-0 ^ 

Z K i D iZ :ci 

i 

e 

\ T n Ax 0 

X K i Di Zj (Ci -Ci) 

Z K.D.zîc? 


\ i ) 


v i y 


( 21 ) 


2.3 Membrane potential 


Introducing eq 13 and eq 21 into eq 1 yields the following expression for the 
membrane potential arising through a porous membrane separating two multi-ionic mixtures 


À'F 

m 


kj , 
— — ln 

z,e 


f . -o V 

rfci 

k B T 

f x — ' Ax 

'Z K.D.z^c, 

i 

— o A 
-Ci) 

ln 

( ,-o \ 

Z K t D iZ fc, 

i 

—Ax n 

l2c, 

e 

ÏK,D lZ ï(ïf 

—o 

-Ci) 



v i y V i y 


( 22 ) 


3 Results and discussion 


Although eq 22 was derived for mixtures with any number of ions, we now focus on temary 
mixtures composed of a symmetric electrolyte (KC1) mixed with an asymmetric one (CaCl 2 ) 
with a common anion. As mentioned in the previous section, we limit ourselves to the 
situation where diffusion drives ail ions towards the same direction. For ail calculations 
presented here, a factor 2 was set for the concentration ratio of each ion, i.e. c° = 2cf x . 


3.1 Influence of the steric, electric and dielectric effects on membrane potential 


Let us define a as the ratio of cation external concentrations 


a = 


(23) 


Figs. 1-4 show the variation of the membrane potential and its components (i.e. the 
différence in Donnan potentials and the diffusion potential) with the concentration of the 
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concentrated compartment for various membrane volume charge densities (X), pore radii (r p ), 
effective dielectric constants inside pores (e p ) and membrane dielectric constants 
respectively. Results shown in Figs. 1-4 were obtained by setting equal cation 
concentrations, i.e. a =1. 

The electroneutrality condition outside the membrane pores reads as follows 


Iz c int =0 


(24) 


The external (i.e. outside the membrane pores) chloride concentration is given 
straightforwardly by introducing eq 23 into eq 24 and we can write 

c mt =(2 + a)c m \ (25) 

cr v ' Ca 2+ 



Figure 1: Membrane potential (-•-), différence between Donnan potentials (— ) and diffusion 

0 3 

potential (-) vs. c Ca 2 + for various volume charge densities (X = -3, -2, -1 eq m' ); a = 1; no 
steric effects; no dielectric effects. 
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Figure 2: Membrane potential (-•-), différence between Donnan potentials (— ) and diffusion 
potential (-) vs. c ° Ca i+ for various pore radii (r p = 1 nm, 3 nm, no steric effects); a = 1; X = - 
1 eq m' 3 ; no dielectric effects. 



Figure 3: Membrane potential (-•-), différence between Donnan potentials (— ) and diffusion 
potential (-) vs. c ^ a 2 + for various effective dielectric constants inside pores (s p = 58, 68, no 
Born effects); a = 1; X = -1 eq m' 3 ; no steric effects; no image charges. 
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Figure 4: Membrane potential (-•-), différence between Donnan potentials (— ) and diffusion 
potential (-) vs. c ° Ca 2 + for various membrane dielectric constants (e m = 3, 6, no image 

O 

charges); r p = 2nm; « = 1; X = -1 eq irr ; no Born effects. 

Let us write 


9î int 


r 

= exp 

v 


-eA'V'p 

k B T 


\ 


J 


(26) 


and 



= A- c 


int 

Ca 


2 + 


( 27 ) 


with Ai= 1, a, 2+a for Ca 2+ , K + and Cl", respectively. 
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Considering eqs 2, 26 and 27, the electroneutrality condition inside the membrane pores (eq 
18) can be rewritten as 


. mt 
Ca 


2 + 


Z z i A ^i 


int 



= -X 


(28) 


Considering KC1 / CaCl 2 mixtures, eq 28 can be rewritten this way 


2 K 


mt 

Ca 


2 + 


(iR int f 


+ a/c 


mt 
K + 



- (2 + a)/c mi =0 
cr 


(29) 


3 2 

According to Cardano’s formula for normal cubic équation as ax + bx + ex + d = 0 , the 
three roots ( jci , xi and * 3 ) are given b y 30 


-3 


- c [ + VÂ 


+= 


-q- 


■ Va" 


+ h 


(30a) 


J-q + y/Â 2 ? — d — VÂ 

x 0 =ûâ — + 0)1 — +h 


(30b) 


x 3 = or \ 


1-q + yfX -q-VÂ 
— + 0)1 — + h 


(30c) 


where 
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ah 3 +bh 2 +ch + cl 


q- 

a 

(31) 

i b 

1 3 a 

(32) 

1 V3 . 

<y= — + — i 

2 2 

(33) 

A 2 V 

A — a + 

27 

(34) 


with p 


2 

3 ah + 2 bh + c 


a 


In the case of KC1 / CaCl 2 mixtures, coefficients a, b , c and d are given by a = 2 K 


int 

Ca 


2 + ■ 


X 


b = a/c 111 ;, , c = — , d = -(2 + a)/c in _ and thus 


K 


c mt 

Ca 


2 + 


cr 


h = ■ 


x + 


6/c 


Ca 


2 + 


(35) 


y A" 


a ' 


/r 


int 

AT 


V y 


P = ~ 


■ + 


X 


/ 


12 


mt 

K r 2 + 

Ca 

v y 


V 2c m * jr mt 


Ca 2+ Ca 2 + 


(36) 


y \ 


er 


/C 


int 
K + 


ax 


q = 


V 


X 


int 
K + 


(2 + a)x m 

cr 


y 


108 


\3 int 


int 

Co 2 * 
v y 


y 


Ca 


2 + 


12 


X2 


int 

Ca 2 ' 
v y 


2/c 


int 

Ca 


2+ 


(37) 
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oc’ 


A=- 

4 


f \3 
K+ 

v y 


\2 


adS- (2 +a)* int 

X ^ cr 


54 ■ & 


\3 int /- 
Ca 2+ g ^nt 


\2 


0? + 


:+ 

V J 


Jnt 

Cd* 


1 

H 

54 


a 2 


y V 

*** 

K* 

v y 


\3 


+— 


X 


f Ÿ r int J nt 

jnt Ca 2+ Q? + 

Cc? + 

v y 


(38) 


Three cases can be considered: 

• A >0: JT’" is a real root and both 31"" and 31"" are imaginary roots. In this case only 
the value of 91| nt is considered. 

• A =0: 9i; nt ,31f , 31"" are ail real roots and 31‘ nt = 91‘ nt . 

• A <0: 9i; nt , 91f 9i; nt are ail real roots. 


3.1.1 Low concentration limit 


For very diluted solutions, it cornes from eqs 37 and 38 



(39) 


limA = 

c) m — >0 


f Y 

1 X 

54 c\X\, 

Ca“ + Ca j 


(40) 


Considering eqs 39 and 40 (q«A) the only one relevant root of eq 29 is 
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9t mt 


fiY /6 


V<V 


M 


- 1/2 


int int 

'-' r , 24 ■V, 2+ 
Ca La y 


(41) 


At low sait concentrations, the co-ions (i.e. chlorine ions) are almost completely excluded 
from the membrane pores. As a resuit, the diffusion potential tends to zero“ as shown in 
Figs. 1-4. The membrane potential is then equal to the Donnan contribution, i.e. 
A^ = AT^ — A 1 ?" • Considering eq 26 the following expression can be derived for the low 

concentration limiting value of the membrane potential 


A¥ = 

m 


k T 


ln 


9T 


<R C 


(42) 


According to eq 42 the limiting value of the membrane potential at low concentration is 
independent of the mixture composition (i.e. a). It is worth mentioning that the non-Donnan 
partitioning coefficient ( je'"' ) is almost independent of the membrane / solution interface (0 or 

Av) at low sait concentrations since (i) both the steric partitioning coefficient ( (/)"' ) and the 

excess solvation energy due to Bom effect ( A VF,- j} orn ) do not dépend on the concentration 


(see eqs 4 and 5) and (ii) the excess solvation energy due to image charges (A W^l iage ) is 

virtually the same at both interfaces because the screening of the ion / image interaction is 
negligible at very low ionic strengths. 25 ' 26 Substituting eq 41 into eq 42 and considering 

c " = 2c" therefore leads to = ln(2) ~ -8 9 mV ’ w hi c h is in agreement with 

numerical results shown in Figs. 1 - 4. It can be noted that this value is identical to that 
obtained with a single sait solution of a 2-1 electrolyte." For a given ratio of the extemal 
concentrations, the low concentration limiting value of the membrane potential in multi-ionic 
mixtures is then ruled by counterions with the highest charge number (z.i). 
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3.1.2 High concentration limit 


The effect of the membrane fixed charge is screened at sufficiently high ionic strength. 

V 

Considering — > o , cq 28 can then be rewritten as 


.A, (9t int ) s ' - 0 


(43) 


At high sait concentrations, the interaction between ions and their images becomes negligible 

'0 'At 

which leads to A Wj image ~ A W i image ~ 0 . Moreover, as previously mentioned both steric 

and Bom effects are independent of the membrane / solution interface (0 or Av). Thus we can 
write (considering once again c° = 2c Ax ) 


—Ax 

cj_ 

-0 
c i 



2 


(44) 


It can be demonstrated that both q and A dépend only on a and k]"' at high sait 
concentrations. Moreover, according to eq 3 w' 1 = xf x = jc i at high sait concentrations. As a 

resuit, the roots of eq 43 dépend only on a and and thus we hâve ° = Aa = . 

According to eq 26, it cornes À v Fp x =À x f / ^. The différence in the Donnan potentials 

(AT^ - A'T®) is then zéro. The membrane potential is then equal to the diffusion potential 
and eq 22 becomes 


A'T 


k B T 

f v i Ax 0 A 

T J K ! D i z i (c i -o-)' 

i 

ln 

( x — ' o 0 A 

Z K i D t z? Ci 

i 

e 

i-» Ax 0 

XA,D,zf(c, - d ) 

v i o Ax 

Z K i D i zfc i 


( i ) 


\ i ) 


(45) 
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Let us write 


— int 

Ck + 


— int 

C C a 2+ 


K 


a 


-9Î 1 


K 


Ca- 


(46) 


Considering eqs 44 and 46, eq 45 can be rewritten as 


A l F. = - 


k B T 


TC 

2 K c^- ! Jr ' ~ K cr°cr 

K r,r 


K, 


2 + cr 


9T 1 
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IC 

^• D of + K r D l! - a ^ L ^' +K cr D cr 


f 


K, 


ln(2) 


(47) 


2 + (x~ 


-9T 1 


K 


V 


Qt- 




Eq 47 shows that the high-concentration limit of the membrane potential dépends on the 
mixture composition (through the cation concentration ratio a) unlike the low-concentration 
limiting value (see eq 42). Eq 47 is also helpful to explain the different behaviors observed in 
Figs. 1-4. When the Donnan exclusion is the only exclusion mechanism (this case is treated 
in Fig. 1), K; =1 and = 1 (because A V E 0 — > Odue to the screening of the membrane fixed 
charge) and eq 47 predicts that the high-concentration limiting value is independent of the 
membrane volume charge density (A) as observed in Fig. 1. The same conclusion can be 
drawn for image charges, i.e. the high-concentration limit of the membrane potential is 
independent of the membrane dielectric constant s m , since A W u ~ 0 . This can be checked 

in Fig. 4. On the other hand, steric effects influence the high-concentration limit of the 
membrane potential (see Fig. 2) since 


K, 


K + 




' Ca 2+ 


( 48 ) 
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Finally, the pore dielectric constant (s p ) dependence of the high-concentration limit of the 
membrane potential observed in Fig. 3 can be explained from eq 47 since Born effects lead to 


K 

K + 

K , 

Ca 2 



Air. )=/(£) 

K + ,Born ' ° x P ' 


(49) 


This last point represents a major différence with the high-concentration limit of the 
membrane potential obtained with single sait solutions since it was shown in a previous work 
that the Born effect does not affect the diffusion potential arising though a membrane 
separating single sait solutions. The high-concentration limit of the membrane potential for 
single CaCl 2 and KC1 solutions can be easily deduced from eq 47 by setting a -» 0 and 
a -» °o , respectively. For example, setting a — » 0 in eq 47 yields 


AT 1 


k B T 


m,CaC ^ 


[ K cr D cr ~ K cr D cr ' 
x 2 K cS D c,r +K cr D cr j 


ln(2) 


(50) 


Eq 50 is équivalent to eq 25 derived in ref 25 and it shows that the influence of the Bom 
effects on the high-concentration limit of the membrane potential disappears in the case of 
single sait solutions. 


3.2 Characterization of membrane pore size and dielectric properties of membrane / 
solution System 


It was shown in a previous work that the effective pore size of membranes can be inferred 
from the high-concentration limit of the membrane potential measured with single sait 
solutions (see e.g. eq 50). 25 Fig. 5 shows the variation of the high-concentration limit of A TV, 
versus the composition of a CaCla / KC1 temary mixture (expressed by the ratio 1 / (1+a)) for 
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various pore radii. The limiting cases 1 / (1+a) = 0 and 1 / (1+a) = 1 correspond to the high- 
concentration limit of the membrane potential of single KC1 and CaCl 2 solutions, 
respectively. It clearly appears from Fig. 5 that the limiting value of A l F„, is much more 
sensitive to the membrane pore size for CaCl 2 than for KC1. As a resuit, it can be concluded 
that asymmetric electrolytes are better candidates than symmetric ones for the détermination 
of membrane pore size from membrane potential measurements. 



AT 

m,CaCl ? 


Figure 5: High-concentration limit of the membrane potential vs. mixture composition 
(c° = 2c") for various pore radii (r p =0.5, 0.7, 1 nm); no Born effect. 


In a recent work, Lanteri et al. suggested to combine membrane potential and rejection rate 
measurements (performed with single sait solutions) so as to détermine the dielectric constant 

o 1 

inside pores (as well as the membrane volume charge density). This procedure was 
successfully applied to a thin-film composite membrane. Here we propose an alternative way 
to assess the effective dielectric constant inside membrane pores ( e p ) with no need to measure 
rejection rates. Indeed, as shown in the above section, the pore dielectric constant ( e p ) may be 
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determined from the high-concentration limit of the membrane potential measured with an 
electrolyte mixture (as illustrated in Fig. 6) provided the pore radius is known (from the 
limiting value of A'F,,, measured with a single electrolyte as discussed above). 



A V F 


m,CaCl 


Figure 6: High-concentration limit of the multi-ionic membrane potential vs. mixture 
composition ( c ° = 2c At ) for various effective dielectric constants inside pores (e p = 58, 68, no 

Born effects); r p = 1 nm. 


4. Conclusion 


The membrane potential arising through charged porous membranes separating two 
aqueous solutions of multi-ionic mixtures at identical hydrostatic pressures but different 
concentrations was investigated through the Nernst-Planck formalism. The Donnan potentials 
that develop at the membrane / solution interfaces were computed by taking into account ion- 
size and dielectric effects. 

In multi-ionic mixtures, it was found that the low-concentration limiting value of the 
membrane potential is ruled by counterions with the highest charge number. The high- 
concentration limit of the membrane potential was found to be dépendent of the mixture 
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composition (a), the pore radius (r p ) and the effective dielectric constant of the solution 
confined inside pores ( e p ). The pore dielectric constant dependence of the high-concentration 
limit of the membrane potential is spécifie to electrolyte mixtures. Indeed it was 
demonstrated in the présent work that this dependence disappears in the case of single-salt 
solutions. The dielectric constant of the membrane material (s m ) does not affect the high- 
concentration limit of the membrane potential because the interaction between ions and their 
images is screened and vanishes at high ionic strength. 

It was shown that membrane potential measured at high sait concentration with single-salt 
solutions and electrolyte mixtures could be used to détermine the pore size and the effective 
dielectric constant inside pores, respectively, with no requirement of additional rejection rate 
measurements as it was proposed in previous works. 
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Conclusion 


La complexité des mécanismes de transfert en nanofiltration (NF) nuit à un 
développement optimal de ce procédé sur le plan industriel. C’est pourquoi les recherches 
visant à relier les caractéristiques d’un nano filtre à ses propriétés de transfert et à mettre au 
point des outils fiables de modélisation du transport, s’intensifient depuis quelques d’années. 

Mon travail a consisté à étendre le modèle de transport SEDE (Steric, Electric and 
Dielectric Exclusion model) à l’étude du transport d’un mélange binaire (2 ions) et ternaire (3 
ions) sous l’action d’un gradient de concentration. L’objectif visé était de développer un outil 
permettant de déterminer les propriétés électrique et diélectrique du système membrane / 
solution à partir d’une grandeur autre que le taux de rejet et caractéristique du transport en 
NF. Cette grandeur est le potentiel de membrane. Celui-ci a tout d’abord été étudié dans le 
cas de solutions binaires (2 ions). Pour une membrane neutre, les résultats obtenus montrent 
que les phénomènes diélectriques (effet diélectrique de Born et forces images) n’ont aucun 
effet sur le potentiel de membrane à la différence de l’effet stérique qui lui, augmentent sa 
valeur. Dans le cas d’une membrane chargée, le potentiel de membrane se situe entre deux 
valeurs limites, celle correspondant au potentiel de Nemst (aux faibles concentrations) et 
celle correspondant au potentiel de diffusion (aux fortes concentrations). A faible 
concentration, les co-ions étant totalement exclus du pore, le potentiel de diffusion s’annule. 
Le potentiel de membrane devient égal à la différence des potentiels de Donnan aux 
interfaces (c’est-à-dire au potentiel de Nemst). A concentration élevée, l’effet Donnan et les 
forces images sont écrantées, les effets stérique et diélectrique de Bom sont indépendants de 
la concentration, il en résulte que les potentiels de Donnan aux interfaces sont identiques. Le 
potentiel de membrane est alors égal au potentiel de diffusion, lequel dépend uniquement des 
coefficients de diffusion libre des ions et des coefficients de gène à la diffusion. Ces derniers 
ne dépendant que du rayon de Stokes des ions et du rayon de pore, la mesure du potentiel de 
membrane à concentration élevée pourrait constituer une méthode alternative pour déterminer 
le rayon de pore de membranes de NF. 

Il est clair que la densité de charge volumique et la constante diélectrique à l’intérieur 
des pores dépendent de l’environnement physico-chimique dans lequel la membrane 
fonctionne. A l’inverse, il est raisonnable de supposer que la constante diélectrique de la 
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membrane est une propriété intrinsèque du matériau membranaire. Ainsi, ni la densité de 
charge volumique, ni la constante diélectrique dans les pores ne peuvent être estimées à partir 
de mesures seules de potentiel de membrane. La détermination de ces deux grandeurs 
nécessite par conséquent d’effectuer des mesures indépendantes telles que des mesures de 
taux de rejet. 

Une partie de mon travail a également consisté à tester la cohérence globale du 
modèle SEDE en analysant deux grandeurs caractéristiques du transport : le potentiel de 
membrane et le taux de rejet ionique. Les expériences ont été menées avec une membrane de 
polyamide et des solutions de chlorure de calcium à différentes concentrations (0.30 à 
llmol.m' ). Le rayon de pore de la membrane et le rapport épaisseur sur porosité ont été 
déterminés respectivement, à partir de mesures de potentiel de membrane relevées à 
concentrations élevées et de perméabilité à l’eau. Un rayon de pore moyen de 0.83 nm et un 
rapport épaisseur sur porosité de 8.24 pm ont été obtenus. Il est à noter que, dans le cas d’un 
mélange binaire (2 ions), le potentiel de membrane ainsi que le taux de rejet d’un sel peuvent 
être décrits de façon satisfaisante par plusieurs couples de valeurs de densité de charge 
volumique (X) et de constante diélectrique à l’intérieur des pores (e p ), parce que l’effet 
Donnan comme l’effet diélectrique de Bom contribuent à l’exclusion d’un électrolyte. 
L’ensemble des couples (X, s p ) reproduisant les mesures de potentiel de membrane fut tout 
d’abord déterminé. Parmi ceux-ci, l’un s’est avéré décrire correctement les taux de rejet 
expérimentaux. Il a été montré que la membrane de polyamide est chargée négativement dans 
la solution de chlorure de calcium et que la densité de charge volumique augmente, en valeur 
absolue, avec la concentration en sel, ce qui suggère une adsorption des ions chlorure à la 
surface des pores. La constante diélectrique dans les pores s’est avérée inférieure à celle de la 
solution externe et indépendante de la concentration en sel. 

La même démarche fut adoptée en considérant seulement les phénomènes d’exclusion 
stérique / électrique, la densité de charge volumique étant cette fois le seul paramètre 
d’ajustement. A la différence du modèle SEDE, la théorie classique s’est révélée incapable de 
décrire à la fois les mesures de potentiel de membrane et de taux de rejet. Ces résultats 
démontrent d’une part, la cohérence globale du modèle SEDE et d’autre part, l’importance de 
phénomènes d’exclusion autres que stérique et électrique, probablement diélectriques, dans le 
cas de membranes de NF. 
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L’étude théorique du phénomène de potentiel de membrane a été étendue à des 
mélanges ternaires (3 ions). Dans le domaine des faibles concentrations, il a été montré que le 
potentiel de membrane tend vers une valeur limite qui ne dépend ni de la composition de la 
solution, ni de la constante diélectrique à l’intérieur des pores et que sa valeur est fixée par 
les contre-ions de plus grande valence. A concentration élevée, le potentiel de membrane tend 
également vers une valeur limite qui dépend cette fois de la composition de la solution. Dans 
ce domaine de concentrations, la charge de surface est écrantée ainsi que l’interaction des 
ions avec leur propre charge image. Le potentiel de membrane est alors indépendant de la 
charge membranaire et de la constante diélectrique du matériau membranaire. En revanche, 
lorsque les effets stérique et diélectrique de Born sont considérés, la limite atteinte par le 
potentiel de membrane dépend de la taille des pores et de la constante diélectrique à 
l’intérieur des pores. La dépendance vis-à-vis de la constante diélectrique à l’intérieur des 
pores est spécifique des mélanges multi-ioniques. En effet, il a été démontré que cette 
dépendance disparaissait dans le cas de solutions binaires. 

Le comportement du potentiel de membrane à concentration élevée est 
particulièrement intéressant puisque sa mesure avec des solutions binaires et ternaires 
permettrait de déterminer la taille des pores et la constante diélectrique à l’intérieur des pores, 
respectivement, sans avoir besoin d’effectuer des mesures supplémentaires de taux de rejet 
d’ions. 
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Perspectives 


Les perspectives à ce travail de recherche s’inscrivent logiquement dans la continuité 
des études présentées dans ce mémoire. 


Suite aux résultats obtenus lors de l’étude théorique du potentiel de membrane en 
présence de mélanges ternaires (préparés à partir de deux sels ayant un ion commun), il 
semble particulièrement intéressant d’effectuer des mesures de potentiel de membrane à 
concentration élevée afin de pouvoir déterminer la constante diélectrique à l’intérieur des 
pores. Ces expériences pourraient être menées avec différents mélanges ternaires, en faisant 
varier la nature des sels et leur proportion afin de pouvoir étudier l’effet du sel sur la 
constante diélectrique de la solution à l’intérieur des pores. 

Il serait également important de comparer les valeurs de la constante diélectrique à l’intérieur 
des pores déduites de la modélisation à des valeurs expérimentales, obtenues par exemple, 
par spectroscopie d’impédance. Cette technique consiste à appliquer, sur une large gamme de 
fréquences, une différence de potentiel alternative à travers une membrane séparant deux 
solutions et à mesurer le courant alternatif résultant. En modélisant le système membrane / 
solution par un circuit équivalent approprié, il devrait être possible de déduire la constante 
diélectrique de la solution à l’intérieur des pores d’une membrane de NF. 

Afin d’estimer au mieux les qualités prédictives du modèle SEDE, il serait également 
nécessaire de déterminer la densité de charge de surface des nanopores. Habituellement, le 
calcul de la densité de charge volumique est effectué à partir de mesures de potentiel 
d’écoulement tangentiel en considérant que la densité de charge électrocinétrique à l’intérieur 
des nanopores de la couche active est identique à celle sur la surface externe. La prise en 
compte, dans le modèle SEDE, du phénomène de régulation de charge n’est pas vraiment 
raisonnable car elle introduirait des paramètres d’ajustement supplémentaires. 

En revanche, le développement de techniques transmembranaires transitoires associé à 
l’analyse théorique des signaux correspondants semble incontournable afin de pouvoir 
accéder aux propriétés de la couche active seule. 
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Annexe 1: Module de 
potentiel de membrane 


La cellule de mesure du potentiel de membrane est représentée sur la figure 16. Les 
électrodes Ag/AgCl sont réalisées par dépôt anodique d’AgCl sur fil d’argent. La membrane 
tubulaire est coupée dans le sens de la longueur, les couches support sont détachées au plus 
de la couche active, et la membrane est rincée plusieurs fois à l’eau distillée (g< 1 |lS. cm ’). La 
membrane est ensuite équilibrée (pendant 12 heures) dans la solution diluée laquelle est 
préparée avec de l’eau de distillée (g< 1 jiS.cm' 1 ). 

Pour finir la membrane est disposée entre deux compartiments remplis d’un même 
électrolyte dans un rapport de concentrations égal à 2 (couche active dirigée vers la solution 
concentrée). Une agitation vigoureuse est réalisée à l’aide d’un agitateur magnétique afin de 
limiter au maximum la polarisation de concentration. Les électrodes sont ensuite plongées 
dans chaque compartiment. Une fois le potentiel de cellule (A'P cell ) mesuré et enregistré, la 

mesure est réalisée une seconde fois en inversant les électrodes afin de s’affranchir de leur 
asymétrie. Le potentiel de cellule sera alors la moyenne arithmétique des deux mesures. 



Figure 16 : Module de mesure du potentiel de membrane 
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Annexe 1 : Module de potentiel de membrane 


Le potentiel de membrane (A^) s’obtient en soustrayant le potentiel de concentration 


k n T 


(A'P^ =— — ln(2) = 17.8mV), traduisant l’influence de la concentration sur le potentiel 


d’électrode, au potentiel de cellule moyen ( A v F cell ) . 


A'f' = A'L 


cell 


A'F 


élect 


eq 71 
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Annexe 2: Module de 
nanofïltration 


Le module de nanofiltration avec lequel nous avons réalisé nos mesures de taux de 
rejet est schématisé sur la figure 17. Une pompe centrifuge fait circuler la solution, contenue 
dans un bac d’alimentation de 25 litres, à travers le module de filtration à un débit volumique 
compris entre 0.4 et 0.5 m .h .Le rétentat est recyclés en continu dans le bac d’alimentation. 
La solution d’alimentation est maintenue à une température constante égale à 25 ± 2°C à 
l’aide d’un échangeur de chaleur. La pression transmembranaire peut être ajustée entre 5 et 
30 bar (± 2 bar) grâce à une vanne disposée en aval du carter. Lors de nos mesures de taux de 
rejet, le perméat est recueilli dans un bêcher, préalablement taré, disposé sur une balance afin 
de déterminer le flux volumique du perméat. 

Avant d’être utilisée, la membrane neuve est disposée dans le carter et conditionnée 
pendant 15 minutes sous une différence de pression transmembranaire de 5 bar en utilisant 
une solution d’alimentation basique (pH 8) puis acide (pH 4) afin de supprimer les résidus 
organiques et minéraux présents dans la membrane. La membrane est ensuite rincée à l’eau 
déminéralisée (a<l|lS.cm ), sous une différence de pression transmembranaire de 5 bar. 

Pour chaque mesure de taux de rejet, un temps de stabilisation d’environ 30 minutes 
est nécessaire. Puis un échantillon de perméat (environ 20 mL) est prélevé, et la valeur du 
flux de perméat enregistré. La conductivité de chaque échantillon de perméat est mesurée afin 
de déterminer le taux de rejet observé en fonction du flux de perméat. 


113 



Annexe 2: Module de nanofiltration 



Figure 17 : Module de nanofiltration tangentielle 
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Résumé 


La complexité des mécanismes de transfert en 
nanofiltration (NF) nuît à un développement 
optimal de ce procédé sur le plan industriel. C’est 
pourquoi les recherches visant à mettre au point des 
outils fiables de modélisation du transport se sont 
s’intensifiées au cours des dernières années. Le 
modèle SEDE (« Steric Electric and Dielectric 
Exclusion ») fait partie de ces outils. Celui-ci 
repose sur des équations de partage aux interfaces 
membrane / solution, incluant des effets stérique, 
électrique et diélectriques, ainsi que sur l’équation 
mésoscopique de Nernst-Planck pour décrire le 
transport dans les pores. L’évaluation de l’aptitude 
prévisionnelle du modèle SEDE passe 

nécessairement par la détermination de ses 
paramètres : taille de pore, densité de charge 
volumique à l’intérieur des pores et constante 
diélectrique dans les pores. Pour déterminer ces 
grandeurs, il m’a été proposé d’étendre 

l’application du modèle SEDE à l’étude du 


potentiel de membrane (différence de potentiel 
électrique s’établissant de part et d’autre d’une 
membrane sous l’action de différences de 
concentration ionique). L’étude théorique du 
potentiel de membrane révèle que la mesure de 
cette grandeur en présence d’un électrolyte binaire à 
concentration élevée, permettrait d’évaluer le rayon 
de pore d’une membrane de NF. Une analyse des 
mesures de potentiel de membrane et de taux de 
rejet ionique effectuées avec une membrane de 
polyamide a permis d’évaluer son rayon de pore, sa 
densité de charge volumique et la constante 
diélectrique à l’intérieur des pores, et de vérifier la 
cohérence globale du modèle SEDE. Enfin, une 
étude théorique du potentiel de membrane avec un 
mélange ternaire (3 ions) révèle que ce type de 
mesures réalisé à concentrations élevées permettrait 
à lui seul d’évaluer la constante diélectrique à 
l’intérieur des pores sans recourir à des mesures 
supplémentaires de taux de rejet. 
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Summary 


The complexity of transfer mechanisms in 
nanofiltration (NF) still harms the full potentialities 
of this process at the industrial scale. That is the 
reason why fundamental studies devoted to the 
development of reliable modeling tools hâve paid 
attention in recent years. The SEDE (« Steric 
Electric and Dielectric Exclusion ») model belongs 
to this category of recently developed modeling 
tools. It is based on partitioning équations that take 
into account the steric, electric and dielectric 
exclusion at the membrane/solution interfaces and 
on the mesoscopic Nernst-Planck équation for the 
description of soluté transport through the 
membrane pores. In order to evaluate the SEDE 
model prédictive capabilities Its parameters, namely 
the pore radius, the volume charge density, the 
dielectric constant of the solution confined inside 
pores and the dielectric constant of the membrane 
material, must be determined. So, we hâve extended 
the application of the SEDE model in estimating 
membrane potential which is defined as the 
electrical potential différence between two 


compartments separated by a membrane and filled 
with two solutions of the same mixture at the same 
température and hydrostatic pressure but different 
concentrations. The theoretical study of the 
membrane potential phenomenon shows that the 
limiting value of the membrane potential obtained 
at high sait concentrations, i.e. the diffusion 
potential, is affected only by steric effects if a 
binary mixture (2 ions) is studied. In other words, 
membrane potential measurements could be used to 
détermine the mean pore size of NF membranes. 
The analysis of membrane potential and ion 
rejection rate data obtained with a polyamide 
membrane has been used to check the global 
cohérence of the SEDE model and it has allowed 
the assessment of the pore radius, the volume 
charge density and the dielectric constant inside 
pores. Finally, in light of another theoretical survey, 
membrane potential in charged membranes 
separating solutions of ternary mixtures (3 ions) 
would allow the assessment of the dielectric 
constant of the solution inside pores. 
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